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ABSTRACT 


This  thesis  investigates  the  feasibility  of  adapting  the  Microflown  PU  match  vector 
sensor  for  underwater  use.  After  testing  the  proper  functioning  of  the  sensor,  the  best 
materials  for  the  capsule  are  detennined  based  on  acoustic  properties.  The  capsule  is 
designed  and  built  by  NPS  staff.  To  predict  the  sensitivity  of  the  encapsulated  sensor,  the 
characteristic  equations  are  modified  to  reflect  the  new  medium.  In  order  to  be  able  to 
predict  the  sensitivity  of  the  encapsulated  sensor,  the  performance  ratio  is  determined 
between  the  sensitivity  in  air  and  in  the  filling  fluid  of  choice.  Temperature  dependency 
is  introduced  in  the  relevant  parameters  to  be  able  to  model  the  sensitivity  at  various 
operating  temperatures.  The  measured  sensitivity  from  the  calibration  report  is  then  used 
to  predict  its  performance  in  castor  oil.  The  theoretical  sensitivity  model  is  verified  by 
experimental  data  gathered  from  calibration  studies  at  NUWC.  The  verified  model  is  then 
used  to  analyze  the  consequences  of  changing  critical  operating  parameters  like  the 
distance  between  the  filaments  and  the  operating  temperature.  Based  on  these 
calculations,  recommendations  are  made  for  a  better  performing  prototype.  Eventually  a 
new  design  is  proposed  that  increases  sensitivity  significantly  and  is  better  adapted  to 
operate  in  the  filling  fluid. 
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I.  INTRODUCTION 


A.  THE  PRINCIPLE  OF  HOT  WIRE  ANEMOMETRY 

Hot  wire  anemometers  have  been  used  in  fluid  mechanics  since  the  late  1800s. 
These  thennal  transducers  are  widely  used  to  study  turbulence  and  unsteady  laminar 
flows.  The  principle  of  hot  wire  anemometry  is  relatively  simple.  A  current  is  passed 
through  a  fine  filament.  When  molecules  flow  by  the  filament  they  will  transfer  heat 
away  from  the  wire.  The  rate  of  heat  transfer  will  be  determined  by  the  flow  speed  and 
flow  temperature,  as  well  as  the  physical  properties  of  the  medium  and  filament.  The 
cooling  of  the  wire  will  change  the  resistance  of  the  filament.  The  variation  in  the 
resistance  can  be  measured. 

There  are  two  basic  modes  of  operation.  The  most  common  is  the  constant  current 
mode,  where  the  current  though  the  filament  is  kept  constant.  The  voltage  difference  to 
accomplish  this  will  vary  with  the  change  in  resistance  of  the  filament.  The  voltage 
difference  thus  becomes  a  measure  for  the  flow  speed  and  temperature.  A  different  mode 
of  operation  is  the  constant  temperature  mode.  This  mode  will  keep  the  filament  at  a 
constant  temperature.  In  order  to  do  this  a  feedback  circuit  will  hold  the  resistance  of  the 
wire  constant  by  varying  the  current  through  the  wire  to  compensate  for  the  heat  transfer 
by  the  flow.  The  current  thus  becomes  a  measure  for  the  flow  speed  and  temperature. 

The  advantage  of  measuring  a  flow  by  using  a  hot  filament  is  the  large  bandwidth. 
A  hot  wire  can  measure  flow  speed  easily  from  DC-flows  to  over  400  kHz.  Furthennore  a 
hot  wire  anemometer  has  very  good  sensitivity  as  well  as  spatial  resolution. 

However,  a  hot  wire  anemometer  is  very  easily  damaged  due  to  the  exposed  thin 
wires.  Additionally,  a  hot  wire  anemometer  is  a  non-linear  system.  This  will  complicate 
its  calibration  process  [1].  Finally,  a  hot  wire  anemometer  needs  a  minimum  flow  speed 
of  about  1  cm/s  to  function  properly. 
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B.  THE  MICROFLOWN  PU  MATCH 

The  Dutch  company  Micro  flown  Technologies  has  taken  a  different  approach  in 
order  to  be  able  to  measure  lower  particle  velocities  using  hot  filaments,  using  two 
parallel  filaments  in  close  proximity.  The  temperature  upstream  is  slightly  less  than  the 
temperature  downstream.  Therefore,  the  downstream  wire  is  heated  more  by  the  upstream 
wire  and  vice  versa  when  the  flow  changes  direction.  The  temperature  difference  between 
the  two  wires  then  becomes  a  measure  for  the  flow  speed  and  direction.  Using  this 
principle  it  can  operate  in  a  flow  range  of  100  nm/s  to  1  m/s  [2],  In  air,  this  means  that 
sound  pressure  levels  (SPL)  from  3.3 dB  to  143 dB  re  20 juPa  can  be  measured.  These  dB 
values  represent  sound  just  over  the  sensitivity  threshold  of  the  human  ear  to  well  over 
the  pain  threshold.  In  seawater  particle  velocities  are  much  smaller  due  to  the  increased 
impedance  of  the  medium. 

The  Micro  flown  Technologies  “PU  match”  sensor  adds  a  miniature  pressure 
microphone  (P)  to  a  single  axis  Microflown  velocity  (U)  sensor  to  fonn  a  package  the 
size  of  a  match.  The  velocity  sensor  measures  the  magnitude  of  the  acoustic  particle 
velocity  along  its  axis.  Combining  the  signal  from  the  velocity  sensor  with  a  pressure 
microphone  enables  the  sensor  to  resolve  whether  the  source  is  in  the  positive  or  the 
negative  direction  along  the  axis  and  to  determine  the  intensity  of  sound  transmitted 
along  the  axis  across  a  wide  bandwidth  (20  Hz-20  kHz). 

C.  POTENTIAL  FOR  UNDERWATER  SOUND  DETECTION 

Although  many  properties  are  the  same  for  both  the  hot  wire  anemometer  and  the 
Microflown  PU  match,  there  are  two  properties  that  motivate  us  only  to  consider  the 
Microflown  PU  match.  In  water,  the  particle  displacement  and  velocity  are  much  smaller 
than  in  air,  but  the  particle  density  is  much  higher.  Since  the  Microflown  PU  match  can 
measure  much  smaller  flow  speeds,  it  is  much  better  suited  for  underwater  sound 
detection  than  traditional  hot  wire  anemometers.  In  addition,  its  beam  pattern  is  a  figure 
eight,  which  makes  unambiguous  directionality  possible  by  adding  a  pressure 
microphone. 
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There  are  several  reasons  why  it  is  interesting  to  examine  the  possibility  to  adapt 
the  Microflown  PU  match  for  underwater  sound  detection.  In  the  first  place,  it  would  be 
able  to  determine  direction  of  sound  without  the  need  for  long  arrays.  This  is  similar  to 
the  advantage  of  traditional  underwater  vector  sensors  based  on  accelerometer 
technology.  Second,  the  response  covers  a  very  large  bandwidth  with  a  single  sensor. 
Moreover  the  Microflown  PU  match  has  very  good  sensitivity  at  very  low  frequencies. 
Since  the  attenuation  of  low  frequency  sound  in  the  ocean  is  much  less  than  that  of  high 
frequencies,  this  capability  could  improve  detection  ranges.  And  finally,  the  small  size  of 
the  sensor  would  enable  integration  in  the  sensor  suite  of  very  small  systems. 

D.  ENCAPSULATION  OF  THE  MICROFLOWN  PU  MATCH 

The  corrosive  properties  of  seawater  make  encapsulation  of  the  sensor  necessary. 
This  provides  challenges  common  to  underwater  sound  transducers.  The  encapsulation 
should  be  some  kind  of  outer  shell  filled  with  a  non-conducting  liquid  (usually  oil  based). 
In  order  to  reduce  reflection  losses  the  impedance  of  both  the  shell  and  the  liquid  must  be 
matched  to  the  impedance  of  seawater.  This  impedance  match  will  also  limit  refraction, 
which  will  introduce  uncertainty  in  the  source  direction. 

It  is  therefore  important  to  carefully  determine  the  materials  used  for  the 
encapsulation.  Since  resources  were  limited,  the  choice  of  outer  shell  materials  was 
primarily  based  on  the  molding  properties  and  less  on  acoustics.  The  choice  of  filling 
fluid  was  however  carefully  based  on  acoustical  properties. 

This  research  is  conducted  to  determine  the  feasibility  to  use  this  sensor  for 
underwater  sound  detection.  At  this  point,  there  is  no  requirement  for  operational 
flexibility.  No  effort  will  be  put  into  enabling  the  sensor  to  operate  in  a  wide  depth  and/or 
temperature  range. 

The  encapsulation  approach  was  designed  and  executed  by  Jay  Adeff  (NPS  staff) 
and  the  details  will  not  be  a  part  of  this  thesis.  A  general  description  will  be  provided. 
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E.  STRUCTURE  OF  THIS  THESIS 

In  this  thesis,  the  perfonnance  of  the  sensor  in  air  at  room  temperature  will  first  be 
evaluated.  This  will  allow  us  to  determine  if  the  sensor  is  working  properly  and  compare 
the  perfonnance  underwater  with  its  perfonnance  in  air,  for  which  it  was  designed.  The 
sensor  performance  underwater  will  then  be  predicted.  The  calibration  tests  to  evaluate 
the  performance  of  the  sensor  underwater  will  be  executed  at  the  Naval  Undersea 
Warfare  Center  (NUWC)  Division  Newport’s  Undersea  Sound  Reference  Division  and, 
due  to  facilities  access,  by  NUWC  staff.  Data  from  the  experiments  will  then  be 
compared  to  the  predictions  and  its  perfonnance  will  be  evaluated.  The  degree  to  which 
alterations  of  the  cunent  sensor  will  influence  its  performance  underwater  will  then  be 
examined.  Finally,  conclusions  will  be  drawn  and  recommendations  will  be  made  for  the 
way  ahead. 
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II.  EVALUATION  OF  THE  MICROFLOWN  PU  MATCH 

IN  AIR 


A.  INTRODUCTION 

The  Microflown  PU  match  is  a  one-dimensional  acoustic  vector  sensor.  It  consists 
of  a  conventional  miniature  Knowles  FG-series  pressure  sensor  and  a  Microflown  Titan 
acoustical  velocity  sensor  element.  The  sensor,  as  shown  in  Figure  1,  is  designed  to 
operate  in  air  between  20  Hz-20  kHz.  The  sensor  is  45  mm  long  and  has  a  width  of 
3.5  mm,  roughly  the  size  of  a  matchstick.  The  upper  sensor  element  is  the  conventional 
pressure  microphone.  The  lower  sensor  is  the  Microflown  Titan  sensor  element.  In  the 
sensor  used  for  this  thesis,  the  Titan  sensor  element  was  oriented  across  the  stem  of  the 
pressure  microphone  element,  resulting  in  a  particle  velocity  axis  at  right  angles  to  the 
stem. 


Figure  1.  Microflown  PU  match  (picture  by  Microflown  Technologies) 

To  verify  its  proper  functioning  and  to  be  able  to  compare  underwater 
performance  to  performance  in  the  medium  for  which  it  was  designed,  the  sensor  was 
first  evaluated  in  air.  The  primary  goal  was  to  establish  the  transfer  function  of  the  sensor 
and  to  see  if  the  pressure  and  acoustical  velocity  sensors  were  coherent.  The  secondary 
goal  was  to  compare  the  performance  of  the  Microflown  PU  match  sub-sensors  to  a  well- 
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known  pressure  sensor;  the  ACO  Pacific  model  7046.  The  beam  pattern  of  the 
Microflown  Titan  sensor  was  also  detennined.  These  measurements  were  all  conducted 
in  the  anechoic  chamber  at  NPS. 


B.  THE  TRANSFER  FUNCTION 

1.  Method 

The  method  used  to  establish  the  transfer  function  in  air  was  based  on  previous 
work  by  Chee  Wee  Ng  and  Jeffrey  Caulk  [3],  [4].  To  determine  the  transfer  function,  the 
sensor  was  placed  in  front  of  an  Austin  AU-15G  sound  source.  See  Appendix  A  for  a 
comprehensive  list  of  all  equipment  used  in  the  measurements.  The  sound  source,  a  guitar 
amplifier,  was  chosen  primarily  for  its  ability  to  produce  loud  noise  over  a  bandwidth  of 
100  Hz-2  kHz.  Other  acoustical  properties  were  of  less  importance  for  this  experiment.  A 
calibrated  ACO  Pacific  pressure  microphone  in  combination  with  an  ACO  pacific  1/2” 
preamplifier  was  installed  above  the  Microflown  at  a  distance  of  3.2  cm.  An  HP  33120A 
Arbitrary  Wavefonn  Generator  was  used  to  produce  white  noise  with  amplitude  of  10VPP. 
The  far  field  condition,  providing  the  range  at  which  wave  front  curvature  may  be 
neglected,  for  a  frequency  of  2  kHz  at  20°  Celsius  [5]  is 


FF  = 


7i  R2 

A 


;r(0.032/n)2 
0. 17m 


(1) 


The  distance  between  source  and  sensor  was  2.04  meter  and  therefore  well  in  the  far 
field. 


Data  acquisition  during  this  experiment  was  conducted  using  the  Mathworks 
Matlab  data  acquisition  toolbox  that  enabled  direct  data  import  from  the  National 
Instruments  cDAQ  (NI  cDAQ).  In  the  system  configuration  used,  the  NI  cDAQ  had  a 
maximum  sampling  rate  of  51.2  kHz.  The  sampling  rates  that  could  be  selected  can  be 
detennined  by  the  following  formula  [6]: 

(2) 

n 
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where  fs  is  the  sampling  frequency,  f„  represents  the  master  clock  rate  (13.1072  MHz) 
and  n  an  integer  between  1  and  31.  Since  the  bandwidth  of  interest  for  this  thesis  is 
20  Hz-2.0  kHz,  a  chosen  sampling  frequency  of  4267  Hz  is  just  above  the  Nyquist 
frequency.  This  corresponds  with/;  =  12  . 


2.  The  Result 


The  estimate  of  the  transfer  function  started  by  establishing  the  frequency  range 
over  which  the  two  sensor  elements  were  coherent.  This  was  done  using  the  Matlab 
command  ms  cohere.  The  coherence  between  pressure  and  velocity  sensor  elements 
increased  significantly  with  higher  source  levels.  The  actual  transfer  function  estimate 
between  the  pressure  element  (p)  and  the  velocity  sensor  (v)  was  calculated  by 


{xp(k)x;xk)) 

\xv(k)x:xk)y 


(3) 


where  Hv{k )  is  the  transfer  function  of  the  kth  frequency  bin  and  Xv (k )  is  the  digital 

Fourier  transform  (DFT)  of  the  time-domain  signals.  For  better  precision  three  datasets 
were  used  to  compute  an  average  transfer  function.  Figure  2  shows  the  average  transfer 
function  for  the  Microflown  PU  match  and  the  coherence  between  the  pressure  element 
and  the  particle  velocity  element.  The  full  analysis  and  computer  code  can  be  found  in 
Appendix  B.  The  output  is  coherent  in  the  bandwidth  80  Hz-2.0  kHz.  Increasing  the 
source  level  had  a  significant  effect  on  the  coherence.  Technical  limitations  of  the 
equipment  used  prevented  further  exploration  of  the  low  end  coherence.  The  transfer 
function  established  during  the  measurements  is  not  smooth.  This  may  be  the  result  of 
reflections  from  the  test  setup  or  some  airflow  in  the  room.  The  amplitude  of  the  transfer 
function  from  the  pressure  sensor  element  to  the  velocity  sensor  element  shows  the  same 
trend  as  follows  from  the  calibration  report.  The  phase  difference  plot  also  follows  the 
prediction  from  the  calibration  report.  From  this  it  can  be  concluded  that  the  Microflown 
PU  match  is  measuring  particle  vibrations  and  performs  as  stated  in  the  calibration  report. 
This  will  allow  us  to  use  the  calibrated  sensitivity  as  a  reference  value  for  future 
calculations. 
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AVERAGED  Coherence  -  Hydroflown:  PRESSURE  to  VELOCITY 
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AVERAGED  Amplitude  of  Transfer  Function  -  Hydroflown:  PRESSURE  to  VELOCITY 
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AVERAGE  Phase  Difference  -  Hydroflown:  PRESSURE  to  VELOCITY 
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Figure  2.  Plot  of  transfer  function  for  the  Micro  flown  PU  Match 


C.  BEAM  PATTERN 

1.  Method 

To  determine  the  beam  pattern  of  the  particle  velocity  element,  the  Micro  flown 
PU  match  was  attached  to  a  Briiel  &  Kjaer  Turntable  system  and  only  the  output  of  the 
particle  velocity  element  was  recorded.  A  Philips  compression  driver  was  fixed  to  a  pole 
at  2.04  m.  This  time  the  ACO  Pacific  reference  microphone  was  not  installed.  A  Stanford 
RS  preamplifier  was  used  for  amplification  and  filtering.  Pre-amplification  was  set  to  lOx 
and  a  1  kHz  High  pass  and  100  kHz  Low  pass  filter  was  used  to  reduce  noise.  An  Agilent 
Waveform  Generator  was  used  to  produce  1000  Hz,  1500  Hz  and  2000  Hz  CW  signals  at 
1.0  Vpp. 

The  turntable  conducted  360  degrees  rotations  with  6  degrees  overshoot  on  both 
sides  to  ensure  a  continuous  rotation  speed  during  a  complete  circle.  Data  was  acquired 
using  National  Instruments  Labview  and  processed  using  Mathworks  Matlab  software. 


2.  The  Result 

During  the  experiment  it  became  apparent  that  the  vibrations  of  the  turntable  had 
a  considerable  effect  on  the  experiment.  Using  a  high  pass  filter  at  1  kHz  and  measuring 
the  beam  pattern  with  frequencies  above  1  kHz  reduced  this  noise  to  acceptable  levels. 
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The  polar  plot  in  Figure  3  shows  the  normalized  response  at  2  kHz  over  a  360°  rotation. 
Although  a  frequency  of  2  kHz  produced  the  clearest  results,  all  three  frequencies  showed 
the  expected  figure  eight  beam  pattern  provided  in  the  data  sheet. 


Beam  Pattern  -  (=2000  Hz 
90 

1 


Figure  3.  The  beam  pattern  of  the  Micro  flown  Titan  in  air 
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III.  SELECTION  OF  THE  MATERIALS  USED  FOR 
ENCAPSULATION 

A.  INTRODUCTION 

The  Microflown  Titan  velocity  sensor  element  uses  two  uncoated  platinum 
filaments  to  measure  the  flow  vector.  Since  the  conducting  filaments  are  uncoated  they 
cannot  be  used  in  direct  contact  with  seawater.  To  use  the  Microflown  PU  match 
underwater  it  therefore  needs  to  be  encapsulated.  This  encapsulation  will  affect  the 
propagation  of  sound  by  reflection  and  refraction.  Since  this  study  is  interested  in  both 
sound  intensity  and  direction,  it  is  important  that  the  capsule  has  as  a  minimal  influence 
on  both  these  properties. 

The  capsule  in  which  the  sensor  will  be  placed  consists  of  two  basic  parts;  the 
outer  shell  and  the  filling  fluid.  For  the  outer  shell  both  the  material  and  the  shape  will 
determine  the  effect  on  sound  transmission  and  direction.  The  filling  fluid  needs  to  be 
non-conducting  and  have  an  impedance  close  to  that  of  average  seawater.  Other  selection 
criteria  were  the  cost  and  availability.  Since  the  goal  of  the  study  is  to  determine  the 
feasibility  to  develop  an  underwater  vector  sensor  based  on  Microflown  technology,  it 
was  considered  acceptable  and  efficient  to  work  with  materials  with  less  than  optimal 
properties. 

The  capsule  was  designed  and  produced  at  the  NPS  by  J.  Adeff.  This  chapter  will 
explain  the  choices  made  in  the  design  and  some  of  the  challenges  during  manufacturing. 

B.  THE  OUTER  SHELL 

The  shell  was  constructed  using  Devcon  Flexane  80  urethane.  This  material  is 
readily  available  and  often  used  for  hydrophones.  The  capsule  has  the  shape  of  a  cylinder 
with  a  hemispherical  end  cap  (see  Figure  4  and  Figure  5).  This  guarantees  a  unifonn  shell 
thickness  to  all  sides.  Furthermore  it  ensured  equidistance  to  the  center  axis  from  all 
sides. 
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Figure  4.  Capsule  (designed  by  Jay  Adeff) 
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Figure  5.  Complete  encapsulated  sensor  along  cm  scale 
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c. 


THE  FILLING  FLUID 


Castor  oil,  Baker  DB  grade,  is  widely  used  in  hydrophones  because  of  its 
excellent  acoustic  properties  and  compatibility  with  other  transducer  materials  [7]. 
Furthermore,  this  castor  oil  is  readily  available  at  limited  cost.  Figure  6  shows  how  the 
characteristic  impedance  of  castor  oil  (Baker  DB  grade)  matches  the  characteristic 
impedance  of  seawater.  Over  the  entire  range  of  0°C-25°C  the  impedance  mismatch 
is  <10%. 


D.  SOME  CRITICAL  ADAPTATIONS 

The  current  Micro  flown  Titan  sensor  element  has  been  designed  for  use  in  air 
only.  The  sensor  uses  platinum  filaments  at  temperatures  around  400°  C  to  optimize 
performance.  The  physical  properties  of  a  liquid  can  be  very  different  from  a  gas.  For 
instance  the  liquid  has  a  boiling  point  at  which  bubbles  will  appear.  This  will  greatly 
decrease  the  sensor  performance  and  probably  reduce  the  life  cycle.  For  castor  oil,  the 
boiling  point  is  313°C  at  atmospheric  pressure.  Due  to  increased  heat  transfer  the  Titan 
sensor  element  will  operate  at  a  slightly  lower  temperature  than  in  air  (the  power 
provided  to  the  sensor  will  remain  the  same),  it  will  still  be  close  to  and  probably  over  the 
boiling  point  of  our  fdling  fluid.  To  prevent  bubbles  from  appearing,  the  operating 
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temperature  needs  to  be  considerably  under  the  boiling  point.  Furthermore  the  platinum 
used  in  the  filaments  operates  as  a  catalyst  in  the  dissociation  of  the  oil.  This  will  not  only 
produce  anomalies  as  the  chemical  composition  of  the  oil  changes,  but  it  will  also  cause  a 
honey-like  residue  on  the  filaments  themselves.  This  can  be  prevented  by  coating  the 
filaments  by  a  thin  layer  of  quartz  [8]. 
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IV.  THE  PERFORMANCE  UNDERWATER 


A.  INTRODUCTION 

Now  that  composition  of  the  prototype  is  known,  a  prediction  of  its  perfonnance 
can  be  made.  In  order  to  do  this  some  assumptions  and  estimates  have  to  be  made.  When 
empirical  data  becomes  available  these  can  be  validated  or  corrected. 

The  first  and  most  important  part  of  predicting  the  perfonnance  of  the  sensor 
underwater  is  to  estimate  how  the  thermodynamic  properties  of  the  filling  fluid  will  affect 
the  heat  distribution.  Particle  displacements  due  to  sound  waves  are  much  smaller  in  a 
liquid  than  in  a  gas.  The  density  of  a  liquid  is  much  greater  than  that  of  a  gas,  so  a  lot 
more  molecules  are  available  to  transfer  the  heat.  Furthermore  the  heat  capacities  of  the 
two  media  are  very  different.  On  the  other  side,  the  diffusion  speed  needs  to  remain  large 
compared  to  the  forced  convection  caused  by  the  sound  wave,  so  that  the  influence  of  the 
particle  displacement  due  to  the  sound  wave  can  be  treated  as  a  small  alteration  to  the 
temperature  profile.  The  start  of  the  prediction  is  therefore  to  calculate  how  all  these 
factors  influence  the  performance. 

The  second  part  of  the  prediction  is  to  estimate  how  the  capsule  will  affect  the 
performance.  What  reflection  losses  are  expected  as  sound  travels  through  the  shell  and 
then  through  the  filling  fluid?  Is  it  possible  to  estimate  the  way  that  sound  direction  will 
be  altered  by  refraction? 

This  chapter  will  first  address  the  influence  of  the  filling  liquid  on  the 
performance  of  the  Micro  flown  Titan  velocity  sensor.  Thereafter  the  influence  of  the 
packaging  will  be  estimated.  The  conclusion  will  be  an  estimate  of  the  sensor 
performance  underwater.  This  will  then  be  compared  to  empirical  data  in  order  to  verify 
theory. 
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B. 


THE  PERFORMANCE  PREDICTION  IN  CASTOR  OIL 


1.  The  Temperature  Difference  between  the  Two  Filaments 

The  sensitivity  of  the  Microflown  sensor  is  based  on  the  temperature  difference 
that  occurs  when  a  flow  passes  two  heated  filaments.  The  temperature  can  be  calculated 
using  the  heat  equation 

pcp(dtT  +  vVT)-kV2T  =  Q,  (4) 

where 
p  =  Density 

c  =  Specific  heat 

T  =  Temperature 
k  =  Thennal  conductivity 
Q  =  Generated  heat 

The  temperature  difference  will  create  a  difference  in  the  resistance  of  the  wires. 
The  change  in  resistance  is  therefore  a  measure  for  the  flow  rate.  A  detailed  calculation  of 
the  temperature  difference  between  the  wires  due  to  sound  waves  can  be  found  in  [9]. 
There,  the  term  vVT  (the  convective  term  in  Equation  (4)  )  is  treated  as  a  perturbation, 
because  in  air  the  diffusion  speed  is  large  compared  to  the  forced  convection.  An 
important  parameter  is  the  heat  diffusion  coefficient  D  =  k  /  pcp.  For  air  at  400°  C  (the 

operating  temperature  of  the  Titan  sensor  element  in  air)  Dair  «  6.5-10  5m2s  1 .  For  castor 
oil  at  3 13°C  (the  maximum  operation  temperature  in  castor  oil)  D  «  3.2  -10  &m2s  1 .  In  the 
current  sensor,  with  a  distance  between  the  wires  (a)  of  100  microns,  the  ratio  between 
the  forced  convection  (v)  and  the  diffusion  speed  (D/a)  at  SPL  « \40dB re  1  pPa 
compares  as 


—  =  3.14 -10  2  «1  and 


D  .  /  a 


^ castor  oil  ^  ^ 


=  2.39-10  2  «1 
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This  shows  that  when  the  sensor  operates  in  castor  oil  the  diffusion  speed  is  still  much 
larger  than  the  forced  convection.  This  implies  that  the  convective  term  can  still  be 
treated  as  a  perturbation  if  the  sensor  is  evaluated  in  oil. 

Including  the  effect  of  the  heat  capacity  of  the  sensors,  the  temperature  difference 
is  given  by  [9] 


AT  =  eilnf 


2  va  P  1 
D  2  7tklv  ifa 


(5) 


where: 


/ 


fo 


f  =  Frequency  [i/z] 


fo  = 


D 


2  na" 


y  ...  2ftfx2 
Jx  D 


fhc  =  ^  ,  j—  =  Second  corner  frequency 


lnLh  (pcX 


h  =  Thickness  of  the  filament 


k  =  Thennal  Conductivity  Coefficient 


W 

m  ■  K 


l  =  Length  of  the  wire  [  m  ] 


v  =  Particle  velocity 


m 

s 


P  =  Power  [w] 
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a  -  Distance  between  wires  [/«] 

L  =  Width  of  the  filament 

K{  (x)  =  Modified  Bessel  function  of  the  second  kind  of  order  1 


The  frequency  /  in  the  first  tenn,  as  stated  in  the  article,  was  replaced  by  the 


dimensionless  frequency  /  ,  assuming  that  it  was  an  error  in  the  article.  This  change  will 
not  impact  the  resulting  analysis  of  ratios  as  this  term  drops  out.  Calculating  the  ratio 


AT.. 


AT. 


oil 


will  eliminate  constant  parameters  from  the  equation,  reducing  uncertainty  with 


regards  to  their  values.  Since  the  sensitivity  is  proportional  to  the  temperature  difference 
between  the  wires,  this  ratio  can  be  used  to  scale  the  sensitivity  function  and  thus  predict 
the  perfonnance.  The  ratio  can  be  broken  down  to 


In  order  to  compare  the  performance  in  air  to  the  performance  in  castor  oil  the 
calculations  will  be  based  on  the  same  amount  of  acoustic  intensity.  The  ratio  between 
particle  speed  amplitudes  (v)  will  therefore  be  based  on  the  same  IL. 


2.  The  Total  Sensitivity  Response 

The  total  frequency  response  of  the  Microflown  can  be  approximated  by  [9] 


|Ar|  =  AT(0) 


1 


1 


(?) 


where 


fd  =  Apparent  comer  frequency  =  some  arbitrary  scaling  factor  *  —  ,  and 

a 
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Ar(o)  = 


va  P 


D  2 nkl. 


(  \ 
.  na 
In  — 
21 

\  V  -v  y 


+  0.577 


This  approximation  will  make  it  much  easier  to  model  the  frequency  response. 
Since  |  A7"|  is  proportional  to  the  actual  response  of  the  Microflown,  the  performance  of 

|A7l 


the  sensor  in  castor  oil  can  be  predicted  by  calculating  the  ratio 


A71 


-.  This  yields 


oil 


f 

1+  \ 

f 2 

1+  ■' 

I^^L'r  _  Ah7  Vair  & oil  \ 
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fd oi, 
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i 
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f 
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(8) 


Using  this  ratio,  the  expected  performance  of  the  encapsulated  Microflown  Titan 
sensor  is  shown  in  Figure  7. 


Frequency  [Hz] 

Figure  7.  The  sensitivity  ratio  in  castor  oil  at  3 13°  C 


Figure  8  shows  that  the  sensitivity  ratio  of  the  Microflown  Titan  element  in  castor 
oil  at  313°C.  At  frequencies  >  1 00  Hz  the  approximation  predicts  less  sensitivity  in 
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castor  oil  than  the  analytic  result  of  the  temperature  difference  ratio.  At  1 000 Hz  the 
sensitivity  ratio  predicted  by  the  approximation  is  twice  the  sensitivity  predicted  by  the 
analytic  result. 

To  be  able  to  predict  the  response,  the  calibrated  sensitivity  in  air  is  needed.  The 
calibration  report  provides  a  graph  of  the  sensitivity  and  a  formula  to  calculate  it.  In  the 
preceding  section  the  sensor  response  in  air  was  confirmed  to  agree  with  this  formula. 
Figure  8  shows  the  calibrated  sensitivity  of  the  Microflown  PU  match  in  air  as  calculated 
from  this  formula  prior  to  encapsulation. 


velocity  sensor  sensitvity  (High  Gain) 


Figure  8.  Sensitivity  in  air  based  on  calibration  report 


Since  the  comparison  of  the  sensitivity  in  air  and  castor  oil  is  calculated  and  the 
sensitivity  in  air  is  known,  we  can  predict  the  sensitivity  in  castor  oil  by: 


Sensitivity castor  oil 


Sensitivity air 
Ratio\AT\ 


(9) 


3.  The  Reflection  Losses 

Reflection  and  refraction  will  occur  due  to  the  impedance  mismatch.  In  the  case 
of  the  encapsulated  Microflown  sensor  in  seawater,  sound  has  to  travel  through  three 
types  of  medium  before  reaching  the  Titan  sensor  element.  Sound  will  first  travel  through 
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the  seawater,  then  encounter  the  Flexane  80  urethane  shell  of  the  capsule.  The  sound  then 
travels  through  the  shell  and  enters  the  castor  oil  to  propagate  to  the  sensor  element.  At 
each  boundary  some  portion  of  the  sound  is  reflected  and  the  direction  of  the  sound  is 
altered.  Although  the  changes  in  apparent  direction  of  sound  are  obviously  important  in  a 
vector  sensor,  this  thesis  will  not  address  this. 

The  impedance  of  seawater  is  well  documented  and  the  impedance  of  castor  oil  is 
also  well  known.  The  impedance  of  Flexane  80  urethane  is  not  known.  Flowever,  the 
shell  has  a  very  small  thickness  (d)  compared  to  the  wavelengths  (X)  on  the  bandwidth  of 

interest.  Even  at  1 000 //z  the  ratio  —  <  0.01 ,  and  can  therefore  be  considered  transparent. 


The  intensity  transmission  coefficient  is  given  by  [10] 


Tj  = 


4  r2rx 


(C  +  h) 


(10) 


where  r  =  pxcx  =  characteristic  acoustic  impedance  of  the  medium. 

The  measurements  at  NUWC  were  conducted  in  seawater  at  20  degrees  Celsius. 
At  that  temperature 


'  Castor  oil 


1.559  ,0“ 

m 


1.446 -106 


Pa-s 

m 


Therefore,  7)  =0.999,  which  means  that  the  intensity  of  the  sound  is  reduced  by  only 

0.1%  before  the  sound  reaches  the  sensor  element.  For  the  comparison  of  the 
experimental  results  with  the  theoretical  model,  therefore,  reflection  losses  can  be 
neglected. 
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4.  The  Predicted  Performance  of  the  Encapsulated  Sensor 

Using  the  calibration  report  and  the  performance  ratios,  the  predicted  performance 
can  be  calculated. 


Predicted  sensitivity,  no  gain 


Figure  9.  The  predicted  performance  of  the  encapsulated  sensor 


Figure  9  displays  the  predicted  performance  of  the  encapsulated  sensor. 
Comparing  the  prediction  to  the  calibration  report  shows  that  the  expected  sensitivity  is 
far  less.  Furthermore  it  shows  that  the  capsule  is  expecting  to  shift  the  peak  sensitivity  to 
a  lower  frequency.  The  sensitivity  of  the  encapsulated  sensor  drops  off  more  quickly  as 
the  frequency  increases. 

At  NUWC,  experiments  will  be  conducted  to  determine  the  actual  sensitivity  of 
the  encapsulated  sensor.  The  results  of  the  experiments  will  be  compared  to  the  predicted 
sensitivity  in  castor  oil  to  verify  the  theory.  If  the  theory  is  confirmed,  it  can  be  used  to 
determine  sensor  adaptations  to  increase  its  performance  underwater. 
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V.  COMPARING  THE  MEASUREMENTS  TO  THE  PREDICTED 

PERFORMANCE 


A.  INTRODUCTION 

At  NUWC,  two  sets  of  measurements  were  done  using  the  encapsulated 
Microflown  PU  match.  These  measurements  were  some  done  a  few  weeks  apart.  During 
the  first  set,  responses  at  two  different  orientation  angles  were  determined,  separated  by 
45  degrees.  First  the  measurements  were  executed  along  the  presumed  axis  of  maximum 
response.  Thereafter,  the  same  measurements  were  executed  at  an  orientation  of  45 
degrees  to  the  assumed  axis  of  maximum  response.  In  the  second  set  of  measurements, 
only  single  orientation  measurements  were  done.  The  single  orientation  measurements 
were  intended  to  correspond  to  the  maximum  response  axis  of  the  particle  velocity 
sensor.  The  aim  of  the  measurements  was  to  determine  the  responses  at  various  Sound 
Pressure  Levels  (SPLs)  and  to  determine  if  the  sensor  still  exhibited  its  directionality 
between  the  two  measurement  angles. 

B.  UNCERTAINTIES 

The  nature  of  the  sensor  causes  some  very  significant  uncertainties  when 
analyzing  its  perfonnance.  First,  the  operating  temperature  is  not  exactly  known.  In  air, 
the  filaments  are  heated  to  approximately  400°  C .  The  larger  heat  conduction  in  castor  oil 
will  reduce  that  temperature  because  the  power  applied  will  remain  the  same.  Since  the 
boiling  point  of  castor  oil  is  313°C  at  standard  atmospheric  pressure,  this  will  be 
assumed  to  be  the  temperature  at  which  the  encapsulated  Microflown  PU  match  will 
operate.  Bubbles  appearing  at  the  filaments  when  the  castor  oil  starts  to  boil  will  further 
influence  the  experimental  results  in  an  unpredictable  way. 

Second,  the  platinum  will  act  as  a  catalyst  in  the  dissociation  of  the  castor  oil  [8]. 
This  process  will  leave  a  honey-like  residue  on  the  wires,  reducing  the  sensitivity  over 
time.  This  process  is  expected  to  influence  the  measurements  within  hours  of  operation. 

Finally,  the  exact  angle  of  refraction  is  not  known,  because  the  impedance  of  the 

Flexane  80  shell  is  not  known.  More  importantly,  due  to  its  small  size,  the  precise 
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orientation  of  the  Titan  sensor  element  within  the  capsule  is  very  hard  to  determine.  Since 
the  Microflown  PU  match  has  a  classic  figure  eight  beam  pattern,  this  will  influence  its 
apparent  sensitivity  when  insonified  at  oblique  angles. 

C.  THE  EXPERIMENTAL  SETUP 


Figure  10.  Experiment  setup  at  NUWC 


The  experiments  to  detennine  the  sensitivity  of  the  encapsulated  Microflown  PU 
match  were  conducted  at  NUWC  by  John  Whitacre  (NUWC  staff).  The  most  important 
part  of  the  experiments  was  to  detennine  the  sensitivity  of  the  sensor  when  a  source  was 
located  along  the  axis  of  maximum  response.  To  determine  if  the  encapsulated  sensor 
retained  its  directionality,  the  sensor  was  also  placed  at  a  45  degree  angle  to  the  initial 
orientation.  During  the  experiments,  the  encapsulated  sensor  was  submerged  in  an 
anechoic  tank  with  a  sound  source  and  a  well-known  A47  SN  16  reference  hydrophone. 
The  bandwidth  between  50-T0007/z was  swept  both  linear  and  logarithmic  at  SPL’s 
between  90-ll0dB  re  \juPa  .  Below  llOdB  re  1  /.iPa  however,  the  travelling  field  was 
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not  created  properly.  During  the  experiment  the  temperature  of  the  tank  was  20°  C .  From 
the  output  voltage  the  sensitivity  was  calculated  in  dB  re  \V  /  juPa  .  These  values  were 
provided  for  further  analysis  and  comparison  to  the  theory. 

The  theory  suggests  that  the  current  encapsulated  sensor  will  be  very  insensitive 
compared  to  its  sensitivity  in  air.  This  means  that  the  expected  output  voltage  is  very  low, 
even  when  high  amplitude  sound  is  directed  at  the  sensor.  The  results  of  the  first  set 
suggested  that  the  sensor  was  indeed  very  insensitive.  From  the  first  set  of  data  it  seemed 
that  the  sensitivity  of  the  sensor  decreased  as  the  SPL  was  increased.  This  could  indicate 
that  the  sensor  was  being  overdriven.  To  further  examine  this  possibility,  a  second,  much 
more  elaborate  set  of  experiments  were  conducted  three  weeks  after  the  first  set.  In  this 
set  the  SPLs  were  reduced  in  3 dB  steps  to  90 dB  re  1  juPa  .  In  total  45  measurements  were 
conducted. 

D.  THE  RESULTS  OF  THE  SENSITIVITY  TEST 

The  first  test  with  a  clear  response  was  obtained  in  the  first  set,  using  a  SPL  of 
140dB,  no  gain.  This  was  the  second  experiment  conducted  with  the  sensor,  so  the  effect 
of  any  residue  on  the  wires  should  be  small.  This  test  will  also  be  used  to  determine  if  the 
encapsulated  sensor  retains  its  directionality. 


Sensitivity  140  dB,  no  gain 
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In  Figure  11,  the  experimental  sensitivity  is  plotted,  together  with  the  theoretical 
prediction  using  the  approximation  and  using  the  ratio  of  analytic  temperature 
differences.  Clearly  the  encapsulated  Microflown  PU  match  is  more  sensitive  to  lower 
frequencies.  The  measurements  show  that  the  apparent  sensitivity  of  the  encapsulated 
sensor  is  highest  around  60 Hz .  It  then  shows  an  unexpected  jumpy  behavior  in 
sensitivity  throughout  the  measured  bandwidth.  On  the  other  hand,  the  position  of  the 
peak  sensitivity  and  the  general  envelope  of  the  sensitivity  seems  to  roughly  correspond 
to  our  model. 

Furthermore,  in  this  case  the  approximation  seems  to  predict  the  sensitivity 
somewhat  better  than  the  direct  calculation  using  the  ratio  of  the  temperature  differences. 
The  difference  between  the  approximation  and  the  exact  solution  is  shown  in  [9].  That 
the  approximation  appears  to  give  a  better  prediction  of  the  experimental  results  can  be 
the  result  of  a  misalignment  of  the  sensor  in  the  capsule.  The  sensor  has  no  calibrated 
direction  markings  on  it.  The  small  size  makes  it  very  hard  to  determine  the  exact  straight 
angle  to  the  sensor  during  the  encapsulation  process.  Because  the  sensor  is  directional, 
any  misalignment  will  result  in  a  reduced  sensitivity  in  a  stationary  experiment.  A  further 
error  in  the  assessed  reflection  losses  caused  by  the  impedance  mismatch  can  further 
influence  the  apparent  sensitivity.  However,  none  of  these  factors  adequately  explain  the 
somewhat  oscillatory  behavior  of  the  sensitivity  and  the  influence  of  electronic  noise  or 
problems  with  the  calibration  process  cannot  be  ruled  out. 
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Sensitivity  120  dB,  no  gain 


Figure  12.  Sensitivity  at  120dB,  no  gain 

In  Figure  12,  the  measured  sensitivity  is  shown  using  SPL  of  120 dB  re  1  /uPa 
without  gain.  The  sensitivity  over  the  bandwidth  changes  in  a  similar  way  as  with  the 
SPL  of  140dB  re  1  juPa  and  again  the  irregular  sensitivity  is  apparent. 

Experiments  conducted  with  SPL  <110<75  re  1  juPa  gave  even  more  irregular 
results.  Sensitivity  seemed  to  randomly  vary  at  different  SPLs.  Additional  experiments 
showed  that  the  traveling  field  was  not  created  well  at  pressure  levels  <11(M5  re 
1  juPa  .Since  the  acoustic  particle  velocity  cannot  be  inferred  from  the  reference 
hydrophone  when  the  contributions  of  the  standing  wave  components  are  unknown,  these 
experiments  could  not  be  used. 

To  better  understand  the  sensitivity  of  the  encapsulated  Microflown  PU  match  in 
comparison  with  the  model,  all  experimental  data  for  SPL>  1  lOdB  re  1  juPa  were 
converted  to  their  values  without  any  gain  and  then  averaged. 
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0.08 


Average  sensitivity,  no  gain 


The  averaged  results  are  shown  in  Figure  13,  together  with  the  predictions  based 
on  the  approximation  and  the  calculated  temperature  difference  ratio.  This  averaging  over 
eight  experiments  at  SPLs  1 10-1 6QdB  re  1  juPa  shows  a  strong  peak  greater  than  the 
expected  value  around  60FIz  followed  by  clear  oscillatory  behavior. 

This  behavior  could  be  caused  by  electronic  noise.  If  60T/zand  subsequent 
harmonics  are  the  reason  for  the  periodicity  in  the  sensitivity,  the  peaks  and  lows  would 
correspond  with  a  harmonic  frequency  of  60 Hz  .  The  initial  peak  is  in  between  56 Hz  and 
64 Hz .  This  suggests  that  it  is  caused  by  electronic  noise.  In  Table  1  the  difference 
between  the  peaks  and  lows  from  the  averaged  experimental  data  and  the 
60Hz  harmonics  is  shown.  The  percentage  difference  seems  to  imply  great  coherence 
between  the  harmonics  and  the  peaks  and  lows.  Fiowever,  the  step  size  during  these 
experiments  was  logarithmic.  This  means  that  the  step  size  started  with  4 Hz  and 
increased  over  the  bandwidth  to  11 0 Hz .  The  peaks/lows  at  higher  frequencies  are 
therefore  determined  less  precisely  than  the  peaks/lows  at  low  frequencies.  To  investigate 
this  variation  further,  the  raw  output  voltage  of  the  measurement  at  1 40dB  and  120 dB  re 
1  juPa  were  compared.  Based  on  the  20 dB  re  1  juPa  difference  in  SPL,  a  factor  of  10  was 
expected  in  the  ratio  of  output  voltages.  Fiowever,  Figure  14  shows  that  the  actual  output 
voltage  ratio  was  only  about  3  and  then  highly  fluctuating  over  the  bandwidth.  The  output 
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voltage  on  the  reference  microphone  showed  the  same  fluctuations.  Inspection  of  the 
traveling  wave  tube  showed  that,  at  the  time  of  the  experiments,  the  projector  might  not 
have  been  functioning  properly.  Therefore,  it  is  justified  to  assume  that,  besides  traveling 
wave  components,  also  some  standing  wave  components  were  present  in  the  tube.  The 
impact  of  the  introduction  of  standing  wave  components  is  unclear  since  the  system’s 
feedback  mechanism  for  generating  traveling  waves  is  unknown.  However,  as  the  system 
sweeps  through  frequency,  the  large  fluctuations  observed  on  the  reference  hydrophone 
strongly  suggest  issues  with  the  desired  system  response.  Unfortunately  time  constraints 
prevented  repeating  the  experiment  in  a  re-calibrated  tube. 
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Figure  14.  Output  voltage  comparison 
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Further  investigation  to  determine  the  cause  of  the  periodicity  in  the  sensitivity  is  needed. 


60  Hz  harmonic 

peak 

low 

%  difference 

60  Hz 

56.7  Hz 

6% 

120  Hz 

135.7  Hz 

12% 

180  Hz 

197.4  Hz 

9% 

300  Hz 

287  Hz 

5% 

360  Hz 

368.4  Hz 

2% 

540  Hz 

535.7  Hz 

1% 

600  Hz 

606.7  Hz 

1% 

Table  1.  Comparison  peaks/lows  to  60FIz  hannonics 

E.  RESULTS  OF  THE  DIRECTIONALITY  TEST 

To  determine  if  the  encapsulated  sensor  retains  directionality,  a  set  of 
measurements  were  conducted  with  the  sensor  at  a  45  degree  angle  relative  to  the  initial 
angle.  Although  the  exact  refraction  angle  is  not  known,  a  45  degree  angle  should  show  a 
significant  drop  in  sensitivity  at  all  frequencies.  To  determine  if  this  was  the  case,  all 
experiments  at  SPL  >11  OdB  re  1  juPa  measured  at  45  degrees  were  averaged  and 
compared  to  the  averaged  sensitivity  of  the  same  SPLs  for  the  initial  orientation.  The 
averaging  was  done  over  four  data  sets  at  both  angles. 
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Figure  15.  Average  percentage  of  sensitivity  loss  at  an  angle  of  45  degrees 

The  results  of  the  comparison  are  shown  in  Figure  15.  Although  the  percentage 
loss  varies  over  the  bandwidth,  it  is  obvious  that  overall  the  sensitivity  is  considerably 
lower  after  a  rotation  of  45  degrees  from  the  presumed  axis  of  maximum  response. 
Therefore,  it  can  be  concluded  that  the  encapsulated  sensor  retains  directivity.  The 
average  loss  was  about  45%.  The  expected  drop  in  sensitivity  would  be  30%  from  the 
measured  beam  pattern  (see  Figure  3).  The  difference  could  be  due  to  additional 
reflection  losses.  An  alternative  cause  is  that  the  initial  alignment  was  not  precisely  along 
the  maximum  response  axis.  If  instead  the  two  measurements  were  at  angles  30  degrees 
and  75  degrees  relative  to  the  maximum  response  axis  (and  45  degrees  relative  to  one 
another),  then  this  drop  in  sensitivity  would  not  be  unexpected.  The  variation  of  the  drop 
in  sensitivity,  however,  is  significant  and  cannot  be  explained  by  a  constant 
misalignment.  This  could  again  indicate  that  the  traveling  wave  was  not  entirely  the  same 
at  the  position  of  the  Micro  flown  during  the  experiments  at  the  two  orientation  angles.  To 
determine  the  actual  beam  pattern,  additional  measurements  are  needed,  where  the  angle 
is  changed  with  small  increments  and  more  data  sets  are  gathered. 
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F. 


EXPERIMENTAL  CONCLUSION 


The  data  shows  that  the  encapsulated  sensor  is  able  to  detect  sound  underwater. 
The  envelope  of  the  sensitivity  roughly  agrees  with  predictions  using  the  calibrated 
sensitivity  in  air  and  the  calculated  sensitivity  ratio  or  the  proportional  ratio  in 
temperature  differences  between  the  filaments.  The  oscillations  in  the  sensitivity  are  not 
predicted  by  the  theory  and  need  additional  measurements  to  investigate.  Calibration 
problems  in  the  traveling  wave  tube  might  explain  some  periodicity  in  the  apparent 
sensitivity  of  the  Micro  flown  Titan  sensor  element.  Electronic  noise  seems  to  have 
influenced  measurements  around  60 Hz .  Further  research  is  needed  to  determine  its 
contribution  to  the  apparent  sensitivity  at  higher  frequencies.  The  encapsulated  sensor 
retains  directionality.  Additional  measurements  are  needed  to  establish  the  exact  beam 
pattern. 


32 


VI.  EFFECTS  OF  SENSOR  MODIFICATIONS 


A.  INTRODUCTION 

The  performance  of  the  Microflown  Titan  is  determined  by  the  difference  in 
temperature  between  the  filaments  as  a  result  of  a  passing  sound  wave.  A  higher 
operating  temperature  will  increase  the  sensitivity.  As  previously  stated,  in  air  the 
operating  temperature  of  the  Microflown  Titan  sensor  element  is  approximately  400°  C . 
The  boiling  point  of  castor  oil  is  3 13° Cat  standard  atmospheric  pressure.  Bubbles 
appearing  when  the  liquid  starts  to  boil  will  greatly  reduce  its  perfonnance.  It  is  therefore 
important  to  operate  the  velocity  sensor  element  at  temperatures  below  the  boiling  point. 
Flow  studies  using  hot  wire  anemometry  by  H.  Eckehnann  have  indicated  that  a 
temperature  of  about  half  the  boiling  point  will  make  the  effect  of  bubbles  appearing 
insignificant  [8].  Since  the  sensor  will  be  used  underwater,  increased  pressure  will 
increase  the  boiling  point  and  therefore  also  increase  the  optimum  operating  temperature. 

The  higher  impedance  of  both  seawater  and  castor  oil  compared  to  air  will  result 
in  a  much  decreased  particle  displacement  and  velocity  caused  by  sound  of  the  same 
Intensity  Level  (IL).  The  amount  of  molecules  available  per  unit  volume  to  carry  heat 
from  one  filament  to  another  however  is  much  larger  in  castor  oil  than  in  air.  These 
parameters  translate  directly  into  the  optimal  distance  apart  at  which  the  filaments  should 
be  placed. 

In  this  chapter,  the  results  of  varying  both  the  operating  temperature  of  the 
filaments  and  the  distance  between  the  wires  will  be  evaluated.  This  can  be  used  to 
determine  these  parameters  in  a  next  prototype  with  the  goal  of  improving  performance. 

B.  THE  METHOD 

At  all  times  the  reference  performance  will  be  the  Microflown  Titan  sensor  as  was 
encapsulated.  This  sensor,  when  used  in  air,  has  an  operating  temperature  of 
approximately  400°  C  and  a  distance  between  the  wires  of  1 00///H  .  These  values  will  not 
change  throughout  this  chapter. 


33 


Several  performance  determining  parameters  are  temperature  dependent.  The  first 
step  in  predicting  the  performance  of  the  sensor  over  a  temperature  range  is  therefore  to 
approximate  the  values  of  these  parameters  over  that  temperature  range.  This  was  done 
by  curve  fitting  data  points  using  Microsoft  Excel  2010  for  the  specific  heat  (cp )  and 

thennal  conductivity  ( k ).  This  approximation  was  necessary  because  only  sparse  data 
was  available.  A  formula  to  calculate  the  density  and  sound  speed  in  castor  oil  was  taken 
from  [7].  To  achieve  different  temperatures,  the  power  (P)  needs  to  become  variable. 
From  the  power  in  the  current  sensor,  the  electric  current  (/0)can  be  calculated. 

Assuming  (70)  to  be  constant  and  using 

P  =  ilR ,  (11) 

and 

R  =  R0a(AT),  (12) 

we  obtain 

P  =  lSHtC,(AT),  (13) 

with 

a  =3.93-10  3 C_1  =  temperature  coefficient, 

70  =  9.9  -10  2  A  =  current  through  the  filament. 

The  resulting  formulas  were  integrated  in  Equations  (5)  and  (7).  To  then  calculate 
the  ratios  of  the  temperature  difference  and  sensitivity  response,  Equations  (6)  and  (8) 
were  altered  to  be  able  to  use  a  constant  value  for  a  in  air  and  a  varying  value  for  a  in 
castor  oil.  In  Figure  16  and  Figure  17  the  specific  heat  and  thennal  conductivity  are 
shown.  Only  two  data  points  were  available  to  determine  the  thermal  conductivity  of 
castor  oil.  The  approximation  was  acceptable,  however,  because  the  total  change  over  the 
entire  temperature  range  was  <  6%  and  consistent  with  similar  oil  types. 
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The  approximated  relationships  will  be  treated  as  valid  over  the  arbitrary 
temperature  range  40°C-400°C ,  even  though  the  boiling  point  of  castor  oil  will  make 
these  approximations  invalid.  We  can  do  this,  because  in  our  comparisons  only  the 
temperatures  below  the  boiling  point  are  used. 
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The  temperature  difference  ratio  now  becomes 


ATair  Kil  Pair  Vair 


AToil  Kir  Kil  Voil 


'  a  (air) 


i -,-x, 
/„,  L 


K, 


a/Z fufiil )  ^0  i\]  if  a  (Oil 


1- 


yjf  i(oii)  P]  yjtf  i(oii)  I  1  —  / 


fhc 


K0  |  y/tfuair )  P()  f  a  (air  ) 


.(14) 


while  the  sensitivity  approximation  ratio  becomes 


\ATL 


1  f 

|  / 2 
'l+  \ 

P oil  Vair  Pair  Kil  \ 

1  fi_on  V 

ffoil 

Dair  Voi,  Poil  Kir 

f /!  i 

1  f2 

1  +  x — 

f  \ 

J  he  air  V 

f 

J  d  air 

(15) 


Using  Equations  (14),  (15)  and  (9)  with  temperature  dependency  incorporated  for 
the  castor  oil  parameters,  the  effect  of  altering  the  operating  temperature  and  the  distance 
between  the  wires  can  be  predicted.  The  experimental  results  have  confirmed  that  the 
sensitivity  in  air  can  be  scaled  to  the  sensitivity  in  castor  oil  at  a  different  operating 
temperature.  This  will  be  used  to  detennine  the  sensitivity  of  the  encapsulated  sensor  at 
operating  temperatures  40°  -3 13° Cat  any  given  distance  between  the  wires. 


C.  THE  EFFECTS  OF  VARIATIONS  IN  OPERATING  TEMPERATURE 
AND  DISTANCE  BETWEEN  THE  WIRES 

The  temperature  dependence  of  the  sensitivity  of  the  current  encapsulated  sensor 
is  shown  in  Figure  18.  The  peak  of  the  current  sensitivity  is  0.04  mV  /  Pa  at  85Hz  and  at 
313°C.  The  sensitivity  deteriorates  with  increasing  frequency  and  with  decreasing 
operating  temperature. 
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Sensitivity  prediction  vs  filament  temperature  for  100  microns 


Frequency  [Hz] 

Figure  18.  Temperature  and  frequency  dependence  of  the  sensitivity  at  100  microns 

Reducing  the  distance  between  the  wires  will  decrease  the  sensitivity  at  very  low 
frequencies.  Flowever,  the  sensitivity  will  be  increased  at  higher  frequencies  [11].  The 
crossover  point  of  the  encapsulated  sensor  is  <  5 Hz ,  so  that  sensitivity  is  in  general 
increased  (see  Figure  19). 


Frequency  [Hz] 


Figure  19.  Sensitivity  at  313°C  and  20qm,  40qm,  60qm,  80  qm  and  lOOqm 


The  distance  between  the  wires  cannot  be  decreased  indefinitely,  however.  The 


sensor  is  based  on  two  filaments  through  which  a  current  flows.  As  the  distance  between 
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the  wires  decreases,  the  induced  magnetic  fields  and  electric  forces  will  start  to  obstruct 
proper  performance.  At  what  distance  this  will  start  to  occur  is  not  part  of  this  thesis. 
However,  the  minimum  distance  will  be  assumed  to  be  >20 /urn  .  Castor  oil  and  quartz 
(possible  material  for  coating  the  platinum  filaments)  are  very  good  insulators  with 
dielectric  constants  of  4.7  and  3.8,  respectively,  at  room  temperature. 

Lowering  the  temperature  will  also  lower  the  sensitivity  over  the  whole  frequency 
range  (see  Figure  18).  Still,  it  is  necessary  to  investigate  the  sensitivity  at  lower 
temperatures  in  order  to  prevent  decreased  performance  by  dissociation  and/or  boiling  of 
the  castor  oil.  The  temperature  at  which  the  oil  is  free  of  boiling  is  dependent  on  the 
pressure  and  therefore  the  operating  depth  of  the  sensor.  Finding  the  exact  temperature  at 
which  dissociation  stops  is  not  a  part  of  this  thesis.  Studies  in  this  field  by  H.  Eckelmann 
[8]  however,  suggest  that  below  half  the  boiling  temperature,  dissociation  becomes 
insignificant.  For  castor  oil  therefore,  an  optimum  operating  temperature  lies  somewhere 
between  100°  -200°C  at  standard  atmospheric  pressure. 

The  combined  effect  of  variation  of  the  operating  temperature  and  the  distance 
between  the  filaments  is  shown  in  Figures  20-22.  The  effects  start  to  significantly  alter 
sensitivity  at  <  40 jum  between  the  filaments. 


Sensitivity  prediction  vs  filament  temperature  for  40  microns 


Frequency  [Hz] 

Figure  20.  Temperature  and  frequency  dependence  of  the  sensitivity  at  40  microns 
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Sensitivity  prediction  vs  filament  temperature  for  30  microns 
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Figure  2 1 .  Temperature  and  frequency  dependence  of  the  sensitivity  at  30  microns 


Sensitivity  prediction  vs  filament  temperature  for  20  microns 


Frequency  [Hz] 

Figure  22.  Temperature  and  frequency  dependence  of  the  sensitivity  at  20  microns 

For  313°C  and  200°  C  the  effect  of  changing  the  distance  between  the  wires  is  shown  in 
Figures  23  and  24,  respectively. 
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Predicted  sensitivity,  no  gain  at  313  deg  C 


Figure  23.  Distance  dependence  of  the  sensitivity  at  313°C 


Predicted  sensitivity,  no  gain  at  200deg  C 


Frequency  [Hz] 


Figure  24.  Distance  dependence  of  the  sensitivity  at  200°  C 


From  Figures  23  and  24  it  can  be  concluded  that  decreasing  the  distance  between 
the  wires  has  a  much  larger  effect  than  changing  the  temperature  in  the  ranges  of  interest. 
Based  on  ratio  of  the  temperature  differences  between  the  wires,  decreasing  the  distance 
between  the  filaments  to  40 jum  will  result  in  a  8dB  gain  in  sensitivity  at  both  313°C  and 
200°C.  Decreasing  the  distance  to  20/wzwill  result  in  a  14r/f?gain  for  both 
temperatures.  A  filling  fluid  with  similar  impedance  and  with  a  much  higher  boiling  point 
might  improve  this  gain  even  further. 
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According  to  the  theory,  decreasing  the  distance  between  the  filaments  will  increase 
sensitivity  rapidly  as  the  distance  between  the  wires  gets  very  small.  This  raises  the 
question  as  to  whether  the  theory  still  holds  at  distances  <30 //m  or  whether  other 
influences  become  dominant.  If  the  theory  does  hold,  it  would  be  useful  to  establish  the 
minimum  distance  apart  that  can  be  achieved. 

The  sensitivity  peaks  just  below  1 00  Hz  in  the  region  where  the  spikes  were 
observed  during  the  experiments  with  the  current  sensor.  If  the  sensor  is  already  sensitive 
in  this  region,  reducing  the  distance  between  the  wires  might  increase  sensitivity  to  the 
point  that  the  sensor  is  overdriven.  In  that  case,  a  design  with  two  sensor  elements  could 
be  considered.  One  sensor  element  optimized  for  frequencies  <150 Hz  and  one  optimized 
for  higher  frequencies.  Since  the  sensitivity  peak  does  not  shift  very  far  as  the  distance 
between  the  wires  decreases  it  is  not  necessary  to  add  more  sensor  elements.  The  size  of 
the  elements  (only  a  few  millimeters)  enables  stacking  them  without  any  significant 
increase  in  overall  sensor  size. 

D.  LIMITS  IN  FILAMENT  PROXIMITY 

Although  theory  suggests  that  a  reduction  of  the  distance  between  the  wires  will 
increase  the  sensitivity  significantly,  technological  restrictions  will  limit  the  extent  to 
which  the  wires  can  be  brought  together.  J.W.  van  Honschoten  has  experimented  in  air 
with  filaments  set30///«  apart  [9].  It  therefore  seems  reasonable  to  assume  that  sensors 
with  a  distance  between  the  filaments  of  about  3  0 // w  can  be  produced. 

E.  SENSOR  MODIFICATION  CONCLUSIONS 

Reducing  the  distance  between  the  wires  increases  the  sensitivity  significantly. 
Reducing  the  operating  temperature  will  have  a  relatively  small  negative  effect  on  the 
sensitivity.  Reducing  the  operating  temperature  in  castor  oil  and  coating  the  filaments  is 
necessary  to  avoid  boiling  and  dissociation  of  the  castor  oil.  Literature  indicates  that  this 
effect  is  negligible  when  the  operating  temperature  is  reduced  to  about  200°  C .  The  effect 
of  reducing  the  distance  between  the  wires  increases  as  the  filaments  are  closer  together. 
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This  encourages  experiments  to  determine  the  validity  of  the  underlying  theory  at  these 
distances.  In  future  designs  of  the  sensor  a  stacking  of  two  elements  could  be  considered 
if  the  sensor  is  found  to  be  too  sensitive  at  low  frequencies. 
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VII.  CONCLUSION  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

It  seems  feasible  to  adapt  the  Microflown  PU  match  for  underwater  use. 
Experiments  seem  to  indicate  that  the  theories  to  predict  its  performance  hold  when  the 
sensor  is  submerged  in  oil.  However,  the  oscillations  of  the  measured  sensitivity  and  the 
high  peak  in  measured  sensitivity  around  60 Hz  are  not  predicted  by  our  model  and  need 
further  investigation.  The  influence  of  the  calibration  error  in  the  experimental  set  up  and 
the  influence  of  electronic  noise  need  to  be  determined. 

Based  on  its  acoustic  and  electric  properties,  castor  oil  is  a  good  choice  as  a  filling 
fluid.  Additional  gain  can  be  achieved  by  finding  a  more  stable  filling  fluid  with  similar 
or  better  acoustic  and  electric  properties,  but  a  much  higher  boiling  point.  Operating  the 
sensor  in  a  liquid  instead  of  a  gas  will  require  some  adaptations  to  the  sensor.  The 
operating  temperature  needs  to  be  well  below  the  boiling  point  of  the  oil.  Coating  the 
filaments  will  be  necessary  to  overcome  catalytic  effects  of  the  platinum  wires.  Further 
strengthening  of  the  wires  may  be  required  to  prevent  damage  when  operating  in  a  much 
denser  medium. 

At  the  current  distance  between  the  filaments,  sensitivity  underwater  is  low. 
Decreasing  the  distance  between  the  filaments  will  increase  the  sensitivity  over  the  whole 
bandwidth  of  interest.  A  gain  of  14-2  IdB  could  be  achieved  if  the  distance  between  the 
wires  is  reduced  to  20-40 jum .  Reducing  the  operating  temperature  to  200°  C  does  not 
significantly  reduce  sensitivity  at  these  distances. 

The  encapsulated  sensor  retains  a  level  of  directivity.  The  exact  beam  pattern 
needs  to  be  detennined  in  future  research. 

B.  RECOMMENDATIONS  FOR  FUTURE  DESIGN 

The  next  prototype  should  have  coated  filaments  to  prevent  the  platinum  from 
acting  as  a  catalyst  in  the  dissociation  of  the  oil.  The  operating  temperature  needs  to  be 
well  below  the  boiling  point  of  the  inert  liquid  of  choice.  Based  on  its  acoustic  properties 
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castor  oil  is  a  good  filling  fluid  for  the  capsule.  In  this  case  an  operating  temperature  of 
200°  C  is  advised.  The  wire  distance  should  be  reduced  to  <40 jum  in  order  to  improve 
sensitivity. 


C.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

Future  research  is  necessary  to  investigate  the  oscillations  in  the  sensitivity.  The 
conducted  experiments  should  be  repeated  in  a  re-calibrated  traveling  wave  tube  to 
determine  if  standing  wave  components  are  the  cause  of  the  periodicity  in  the  sensitivity. 
Additional  experiments  are  necessary  to  determine  the  influence  of  electronic  noise. 
Subsequent  research  is  necessary  to  determine  if  the  beam  pattern  of  the  encapsulated 
sensor  maintains  its  proper  response.  Finally,  it  needs  to  be  determined  if  the  current 
prediction  model  will  hold  at  very  small  distances  between  the  wires. 
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APPENDIX  A.  LIST  OF  EQUIPMENT 


Equipment  and  settings  used  during  all  measurements 

Microflown  PU  match  consisting  of: 

-  Microflown  Titan  acoustical  particle 

velocity  sensor 

-  Knowles  FG  series  miniature  sound 

pressure  transducer 

Frequency  range: 

-  Velocity  sensor:  0.1  Hz-20  kHz 

-  Pressure  sensor:  20  Hz-20  kHz 

Microflown  Technologies  MFSC-2 

2  channel  signal  conditioner 

-  Gain  set  to  ‘High’ 

-  Correction  mode  set  to  ‘Off 

Post-processing: 

Mathworks  Matlab  Version  2010a 

Additional  equipment  and  settings  used  to  determine  beam  pattern 

Agilent  type  33220A 

Arbitrary  Waveform  Generator 

-  1000  Hz,  1500  Hz,  2000  Hz  CW 

-  Vppsetto  ‘1.0  V’ 

Briiel  &  Kjaer  Turntable  system  type  9640 

Stanford  RS  preamplifier  model  SR560 

Pre-amplification:  lOx 

High  pass  filter:  1  kHz 

Low  pass  filter:  100  kHz 

Philips  compression  driver 

Data  acquisition  software: 

National  Instruments  Labview  2011 

Sampling  frequency: 

25  Hz 

Additional  equipment  and  software  used  to  detennine  transfer  function 

National  Instruments  Compact  DAQ  USB 

chassis  model  NI  9172 

National  Instruments  Sound  and  Vibration 

DAQ  module  model  Nl  9234  (2  total) 

45 


Austin  AU-15G  Amplifier 

-  Overdrive  set  to  ‘ON’ 

-  Volume  set  to  ‘MAX’ 

Hewlett  Packard  Function  Generator  model 

33120A 

-  Mode:  White  Noise 

-Amplitude:  10VPP 

AGO  Pacific  Calibrated  Pressure 

microphone  model  7046 

ACO  Pacific  V”  pre-amplifier  model  4012 

Data  acquisition  software: 

-  Mathworks  Matlab  2011 

-  Data  acquisition  toolbox 

Additional  equipment  and  software  analyze  data  NUWC 

Microsoft  Excel  2010 
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APPENDIX  B.  EVALUATION  OF  THE  MICROFLOWN  IN  AIR 


A.  THE  FULL  EVALUATION 

The  proper  functioning  of  the  Microflown  PU  Match  sensor  and  all  its  sensor 
elements  was  confirmed.  To  be  able  to  do  this,  first  the  transfer  function  was  established 
between  the  pressure  element  and  the  velocity  element.  Also  the  transfer  function  with  an 
independent  and  well  known  ACO  Pacific  model  7046  microphone  was  established  and 
compared  to  the  initial  results.  Finally  the  Noise  Spectrum  Level  (NSL)  was  determined. 
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B. 


THE  MATLAB  CODE  FOR  THE  TRANSFER  FUNCTION 


%%  Hydroflown  Transfer  Function 

%  Program  to  read  in  Hydroflown  data  and  compute  transfer  function  between 
%  pressure  sensor  and  velocity  sensors 


%%  Clean  up 
clear  all 
cl  c 

close  all 

%%  Setting  up  F  F T 
maxchan  =  6; 
numpass  =  250; 

N  =  1024; 

M  =  numpass  *  N; 

Fs  =  4267; 

NF  FT  =  2A10; 

overlap  =  NF  FT/  2; 

wi  ndow  =  hanni  ng(  NFFT) ; 


%ma ximum  number  of  data  channels  captured 
%number  of  loop  iterations  in  vi  sampling  routine 
%number  of  samples  per  channel  per  pass 
%t  o  t  a  I  sample  length  per  channel 
%s  a  mp I  i  n  g  f  r  e  q 


%%  Channel 
p  r  s  1  =  1 ; 
bl  ul  =  2; 
g  r  n  1  =  3 ; 
p r  s  R  =  4; 
p  r  s  2  =  5 ; 
bl  u 2  =  6; 


Desi  gnat i on 
%  sensor  324  pressure 
%  sensor  324  blu  veloci 
%  unused 
%  A  C  0  pressure 
%  unused 
%  unused 


ty 


Srefl  =  prsl;  %  hydroflown  Pressure  Channel 
S  r  e  f  2  =  bl  ul;  %  hydroflown  Velocity  Channel 


pi  ot  sensor  ='  Mi  crofl  own'  ; 
pi  ot  ref  1='  PRESSURE'  ; 
pi  ot  ref  2='  VELOCI  TY'  ; 
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%%  Loadi  ng  Datasets 


bindata  =  zeros(M,  maxchan);  ^initialize  matrix  to  allocate  memory 


datasets  =  3;  %n u mb e r  of  data  sets  for  averaging 


for  ds  =1:  dat  as  et  s  %c  y  c  I  e  through  number  of  data  sets 


if  ds  ==  1 
bindata  = 

i  mport  data( '  C:  \  Users  \  Marni  x  \  Documents  \  MAT  L  A  B  \  Thesi  s  \  measurement  s  \  Mar20_  WN  _10Vpp_max 
_l_60s,mat’); 

else  if  ds  ==  2 

bindata  = 

i  mportdataf'  C:\Users\Marni  x\Documents  \  MAT  L  A  B\  Thesi  s  \  measurement  s  \  Mar20_WN_10Vpp_max 
_2_60s.  mat  ’  ) ; 


else  if  ds  ==  3 


bindata  = 

i  mport  data{ '  C:  \  Users  \  Marni  x  \  Documentst  MAT  L  A  B  \  Thesi  s\measurements\Mar20_  WN  _10Vpp_max 
_3_60s.  mat1  ) ; 


end 


end 


end 


bindata  =  bindata’;  %c h a n g e  to  row  data 


%%  C  o  mp  u  t  a  t  i  o  n  for  Transfer  Function 
channel 

g  a  mma  _  p  2  v  (  :  ,  d  s )  = 

mscoherefbi  ndatafSrefl,  ;  ),  bi  ndata(Sref2, 

[  P 2 V( ;  ,  ds ) ,  F ]  =  c ps d ( bi  nda t a ( Sr ef  1 , : 

[V(:,ds),F]  =  pwe  I  c  h  (  bi  nda  t  a  ( S  r  ef  2 ,  : 

[P(;,ds),F]  =  pwe  I  c  h  (  bi  nda  t  a  ( S  r  ef  1 ,  : 

Hp  2  v  (  ;  ,  ds)  =  P  2  V  (  :  ,  ds )  .  /  V(  :  ,  ds ) ;  "/ 


for  Sound  Probe  Pressure  to  velocity 

) ,  wi  ndow,  over  I  ap,  NFFT,  Fs) ; 

),  bi  ndata(Sref2,  ;  ),  wi  ndow,  overl  ap,  NFFT,  Fs) ; 
) ,  wi  ndow,  over  I  ap,  NFFT,  Fs)  ; 

) ,  wi  ndow,  over  I  ap,  NFFT,  Fs)  ; 

,  PRESSURE  to  Vel  oci  ty 
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%%  Computation  for  ACO  Reference  to  Pressure  element  channel 


gamma_pACO(  :  ,  ds)  = 

ms  coher  e(  bi  n  d  a  t  a  (  prsR,  :  ) ,  bi  n  d  a  t  a  ( Sref  1,  :  ) ,  wi  ndow,  overl  ap,  NFFT,  Fs) ; 
[  P2ACO(  :  ,  ds) ,  F]  = 

cpsdjbi  ndatafprsR,  :  ),  bi  ndatafSrefl,  :  ),  wi  ndow,  overl  ap,  NFFT,  Fs); 

[P(:,ds),F]  =  pwel  ch(  bi  ndataf  Sref  1,  :  ),  wi  ndow,  overl  ap,  NFFT,  Fs) ; 

F 2 x ( ; , ds)  =  P 2 A C 0 ( ; , ds ) . / P( : , ds ) ;  %  ACO  t  o  PRESSURE 


end 


%%  PLOT  Transfer  functions 

%  VELOCITY  element  to  PRESSURE  element  Transfer  Function 
f  i  gu  r  e(  1) 
subpl  ot  ( 3,  1,  1) 
hoi  d  on 

pi  ot  ( F,  gamma_p2v(  ;  ,  1) ,  '  b'  ) 
yl  abel  ( '  Coherence'  ) 

t  i  1 1  e(  [ '  Coher  ence  •  .plotsensor,  .plotrefl,  to 
,  pi  o  t  r  e  f  2  j ,  1  Font  Wei  g  h  t ' , '  bold’) 

a x i  s  ( [  0  2  0  0  0  0  1.1  ]) 

subpl  ot  ( 3,  1,  2) 
hoi  d  on 

pi  ot(F,  abs(  Hp2v(:  ,  1)),  ’  b'  ) 

yl  abel  ( '  Tr  ansf  er  Function  Amplitude') 

t  i  1 1  e(  [ '  Ampl  i  t  ude  of  Transfer  Function  ■  '.plotsensor,  .plotrefl, 
,  pi  o  t  r  e  f  2  ] ,  'Font  We  ight',  'bold') 

a x i  s ( [  0  2  0  0  0  0  2  ] ) 


subpl  ot  ( 3,  1,  3) 
hoi  d  on 

pi  o  t  ( F ,  rad2deg(  angl  e  ( H  p  2  v  (  ;  ,  1))),  '  b'  ) 

t  i  1 1  e  ( [ '  P  h  a  s  e  Difference  -  .plotsensor,  .plotrefl,'  to 
,  pi  o  t  r  e  f  2  ] ,  '  Font  We  ight',  'bold') 

yl  abel  (' Phase  Difference  (deg)’) 


t  o 
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xl  abel ( '  Frequency  ( Hz) ’ ) 

a x i  s ( [  0  2  0  0  0  -  1  8  0  1  8  0  ] ) 

%  PLOT  averaged  VELOCITY  element  to  PRESSURE  element  Function 
f i  gure( 2) 
s  ubpl  ot  ( 3, 1,  1) 
hoi  d  on 

gamma_p2v=sum(  gamma_p2v,  2)  /  3; 

pi  ot  ( F,  gamma _ p 2 v ,  '  b' ) 

yl  abel  ( '  Coherence' ,  '  FontSi  z e '  ,  14) 

t  i  1 1  e(  ['  AVERAGED  Coherence  -  ,  pi  otsensor,  .plotrefl,'  to 

,  pi  o  t  r  e  f  2  ]  ,  '  FontSi  z  e '  ,  18,  '  Font  Wei  g  h  t  ’  ,  '  bol  d  ’  ) 

a x i  s  ( [  0  2  0  0  0  0  1.1  ]) 

subpl  ot  (3,  1,  2) 
hoi  d  on 

H p 2 v  =  sum]  Hp2v,  2)/3; 
pi  ot(  F,  abs( Hp2v) ,  '  b'  ) 
yl  abel  ( '  A  mp  I  i  t  u  d  e ' ,  '  FontSi  z  e ' ,  14) 

t  i  1 1  e  ( [ '  AVE  RAGE  D  Amplitude  of  Transfer  Function  -  '  ,  pi  otsensor,  , plotrefl, 
to  '  ,  pi  otref2] ,  '  FontSi  ze'  ,  18,  '  Font  We  ight'.'bold’) 

a x i  s ( [  0  2  0  0  0  0  2  ] ) 

subpl  ot  (  3,  1,  3) 

pi  o  t  ( F ,  rad2deg(  angl  e  (  H  p  2  v ) ) ,  '  b’  ) 

t  i  1 1  e  ( [ '  AVE  RAGE  D  Phase  Difference  -  .plotsensor,  .plotrefl,'  to 

,  pi  o  t  r  e  f  2  ] ,  '  FontSi  ze'  ,  18,  ’  Font  Wei  g  h  t  ’  ,  '  bol  d ' ) 

yl  abel  ('  Phase  Difference  ( deg) ' , '  Font  Si  z  e  ’ ,  14) 

xl  abel  ( '  Frequency  ( Hz)  ’ ,  ’  FontSi  ze’ ,  14'  ) 

a x i  s ( [  0  2  0  0  0  -  1  8  0  1  8  0  ] ) 

%  PLOT  PRESSURE  element  to  ACO  microphone  Transfer  Function 
f  i  gu r  e( 3) 
subpl  ot  ( 3,  1,  1) 
hoi  d  on 
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p I  o t  ( F ,  gamma_pACO( :  ,  1) ,  '  b'  ) 
yl  abel { ' Coherence' ) 

t  i  1 1  e(  [ '  Coher  enc  e  -  ,  pi  otsensor,  '  :  ACO  Microphone  to  ‘ 

plotrefl], 'Font  We  ight', 'bold') 

a x i  s ( [  0  2  0  0  0  0  1.1  ]) 

subpl  ot(3,  1,  2) 
hoi  d  on 

pi  ot{  F,  abs(  F2x( : ,  1) ), '  b* ) 

yl  abel  ( '  Tr  ansf  er  Function  Amplitude') 

t  i  1 1  e(  [ '  Ampl  i  t  ude  of  Transfer  Function  -  '  ,  pi  otsensor,  :  ACO  Microphone  to 
plotrefl],  'Font  We  ight', 'bold') 

a x i  s ( [  0  2  0  0  0  0  2  ] ) 

subpl  ot  (  3,  1,  3) 
hoi  d  on 

pi  o  t ( F ,  r  ad2deg( angl  e ( F  2  x ( : ,  1))), '  b’ ) 

t  i  1 1  e  ( [ 1  P  h  a  s  e  Difference  -  , plotsensor,  ACO  Microphone  to  ' 

plotrefl],  Font  We  ight1, 'bold') 

yl  abel  ('  Phase  Di  f  f  er  ence  (deg)') 

xl  abel  ( '  Frequency  ( Hz) ' ) 

a x i s ( [ 0  2  0  0  0  ■  1  8  0  1  8  0  ] ) 

%  PLOT  averaged  PRESSURE  element  to  ACO  Transfer  Function 
f  i  gu  r  e( 4) 
subpl  ot  ( 3,  1,  1) 
hoi  d  on 

gamma_pACO=sum(  gamma_pACO,  2)  /  3; 
pi  ot  ( F,  gamma_pACO,  '  b'  ) 
yl  abel  ( ' Coherence' ) 

t  i  1 1  e(  [' AVERAGED  Coherence  -  , plotsensor,  ACO  Microphone  to  1 

plotrefl], 'Font  We  ight', 'bold') 

a x i  s  ( [  0  2  0  0  0  0  1.1  ]) 

subpl  ot  ( 3,  1,  2) 
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hoi  d  on 


F  2  x  =  s  u  m(  F  2  x ,  2 )  /  3 ; 
plot(F,abs(F2x),'b  ) 
ylabelf'Transfer  Function  Amplitude’) 

t  i  1 1  e(  [’ AVERAGED  Amplitude  of  Transfer  Function  ■  pi  ot  sensor ,  ACO 
Microphone  to  ‘  ,  pi  otrefl],' Font  Weight’ bold') 

a x i  s ( [  0  2  0  0  0  0  2  ] ) 

subpl  ot  ( 3,  1,  3) 
hoi  d  on 

pi  ot  ( F,  r  a d 2 d e g ( angl  e(  F 2 x ) ) ,  ’  b' ) 

t  i  1 1  e  ( [ '  AVE  RAGE  D  Phase  Difference  -  .plotsensor,  ACO  Microphone  to  ' 

,  pi  otreflj,’  Font  Wei  g  h  t ' , '  bold’) 

yl  abel  ('  Phase  Di  f  f  er  ence  (deg)') 

xl  abel  ( '  Frequency  ( Hz) ’ ) 

a x i  s ( [  0  2  0  0  0  -  1  8  0  1  8  0  ]  ) 

%  PLOT  VELOCITY  element  to  ACO  microphone  transfer  function 
f  i  gu  r  e( 5) 
subpl  ot  ( 2,  1,  1) 

F  H  =  H  p  2  v .  *  F  2x;  %transfer  function  to  correct  velocity  to  ACO  pressure  direct 

pi  ot  ( F,  a b s ( FH) ,  ’  b' ) 

yl  abel  {'Transfer  Function  Amplitude') 

t  i  1 1  e(  [ 1  Ampl  i  t  ude  of  Transfer  Function  •  ' ,  pi  ot  sens  or ,  ’  :  ACO  Microphone  to 
' ,  pi  otref2], '  Font  We  ight'.'bold') 

a x i  s ( [  0  2  0  0  0  0  3  ] ) 

s  ubpl  ot  ( 2,  1, 2) 

pi  ot  ( F,  r  a d 2 d e g ( angl  e(  FH) ) ,  '  b' ) 
yl  abel  ('  Phase  Di  f  f  er  ence  ( deg) ' ) 
xl  abel ( '  Frequency  ( Hz) ’ ) 

t  i  1 1  e(  [ '  Phase  difference  of  Transfer  Function  ■  ,  pi  otsensor,  '  ACO  Microphone 

to  ’  ,  pi  o  t  r  e  f  2  ] ,  ’  F  o  n  t  We  i  g  h  t  ’  ,  ’  b  o  I  d  ) 

a x i  s ( [  0  2  0  0  0  -  1  8  0  1  8  0  ] ) 

%%  Determine  Noise  Spectrum  Level 
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bindata  =  zeros(  M,  maxchan) ;  %i  n  i  t  i  a  I  i  z  e  matrix  to  allocate  memory 


bi  ndat  a  = 

i  mpo r  t  d a t  a ( '  C:  \  Users  \  Marni  x  \  Documents  \  MAT  L  A  B  \  Thesi  s\measurements\Mar20_  WN  _10Vpp_max 
_3_60s.  mat'  );  %i  n  i  t  i  a  I  i  z  e  ma  t  r  i  x  to  allocate  me  mo  r  y 

m=0 .  0  5  4  3  3;  %aco  sensitivity 

bindata  =  bindata';  %c  h  a  n  g  e  to  row  data 

noi  se_si  gnal  =  bindata(Sref2,:); 

noi  se_si  gnal  =  noi  se_si  gnal  -  mean(  noi  se_si  gnal  ) ; 

pent  =  p  we  I  ch(  noi  s  e  _  s  i  gnal  ,  wi  ndow,  overl  ap,  NFFT,  Fs) ; 

pe  nf  =pe  nf .  *  ( F  H.  *  c  o  n  j  ( F  H)  )  /  mA2*  1 E  1  2/  4  0  0;  %pressure  equivalent  noise  floor 
computation  re  2 0 u P a 

pent  =  10* I og  10( penf ) ; 

ACO  =  bi  ndataf  prsR,  :  )  ; 

ACO  =  ACO  -  mea  n(  ACO) ; 

[PACO,  F]  =  pwel  ch(ACO,  wi  ndow,  overl  ap,  NFFT,  Fs)  ; 

PACO  =  PACO/mA2;  %  ACO  corrected  to  Pa 
PACO  =  PACO*  1 E  1  2/  400 ;  %  ACO  signal  in  2 0 u P a 

PACO  =  10*1 o g 1 0 (  PACO) ; 
di  f  f  =P ACO-  pent ; 


%%  Plot  Noise  Sped  r  urn  Level  s 

%  PLOT  Noise  Spectrum  Level  measured  by  ACO  microphone  and  PRESSURE  element 
f  i  gur  e( 6) 
s  ubpl  ot  ( 3,  1,  1) 
pi  ot(  F,  pent, '  b' ) 

xl  abel  ( '  Frequency  ( Hz) ' ) 
yl  abel  ( '  NSL  re  20uPaA2/  Hz’  ) 

t  i  1 1  e(  ['  Pressure  Equivalent  Noise  Spectrum  Level  of  .plotsensor, 
plotref2], ’Font  We  ight’,' bold’) 

a x i s ( [ 0  2  0  0  0  0  5  5  ]) 
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s  ubpl  ot  ( 3,  1,  2) 
pi  ot  ( F,  PACO,  '  r ' ) 
xl  abel  ( '  Frequency  ( Hz) ' ) 
yl  abel  ( '  NSL  re  20uPa"2 /  Hz’ ) 

ti  tl  e('  Pressure  Equivalent  Noise  Sped  r  um  Level  of  ACO 
Mi  crophone' ,  1  Font  Wei  g  h  t  ’  ,  ’  bol  cT  ) 

a x i s ( [ 0  2  0  0  0  0  5  5  ] ) 

s  ubpl  ot  ( 3,  1,  3) 
pi  ot ( F,  di  f f ,  b'  ) 
hoi  d  on 

%  p  I  o  t  ( F ,  p  e  n  f ,  '  r '  ) 
xl  abel  ( '  Frequency  ( Hz) ' ) 
yl  abel  ( '  NSL  re  20uPaA2/  Hz'  ) 

ti  tl  e('  Di  fference  in  Pressure  Equivalent  NSL  of  ACO  Microphone  and 
Velocity', 'Font  We  ight'.'bold') 

a x i s ( [ 0  2  0  0  0  -  1  1  0  ] ) 

MATLAB  CODE  FOR  DETERMINING  THE  BEAM  PATTERN 

%%  Cal  c  u  1 1  e  beampattern  on  Microflown  particle  velocity  sensor 
%  frequency  1500  Hz 
%  S  a  mp  I  i  n  g  period  4 .  0  0  e  -  2  s 
cl  ear; 
cl  c 


c 1 c  =x I  sread]  2000_cropped_cl  ean, xl 

tt  =4.  OOe- 2*1 engt  h( clc)  ; 
c 1 c ( :  ,  1)  =clc(  :  ,  1)  -  c 1 c ( 1,  1) ; 
a  mp  =c 1 c  (  :  ,  2 ) ; 
dbamp=20.  *1  ogl0(  amp)  ; 
dbamp=dbamp-  mi  n(  dbamp)  ; 
dbamp=dbamp  /  max]  dbamp) ; 
angamp]  ;  ,  l)=clc(:  ,  1 )  *  ( (  2  *  p  i  )/tt)  ; 
angamp]  ;  ,  2)  =d  b  a  mp ; 


sx'  ); 

%  t  o  t  a  I  t i  me  [  s  ] 

%  st  art  at  t  =0  s 
%  define  a  mp  [  V  ] 

%  convert  to  decibels 
%  s  e  t  mi  n  i  mu  m  t  o  0  d  B 

%  convert  time  to  radians 
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pol  ar(  angamp(  :  ,  1) ,  angamp(  :  ,  2) ) ;  %  plot  beampattern 

titlef'Beam  Pattern  ■  f  =2  0  0  0  Hz'); 
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APPENDIX  C.  MATLAB  CODE  TO  DETERMINE  ACOUSTICAL 
PROPERTIES  OF  CASTOR  OIL 


%%  R h o _ 0 ,  c  Calculation  of  Castor  oil 

%  Intended  to  gain  insight  in  the  acoustic  properties  of  castor  oil  and  seawater  as 
a 

%  function  of  temperture 
%%  Clean  up 
clear 
cl  c 

close  all 
%%  p  a  r  a  me  t  e  r  s 

t  =0 :  0 .  1 : 2  5;  %Temperture  range  of  the  oil  indegCelcius 

p  =2  0.265;  %  Pressure  in  dBar 

SA=35*ones(  1,  I  ength(t) ) ;  %  Salinity  [g/kg]  (constant  in  this  case) 

%%  Calculation  density  and  speed  of  sound  Castor  oil 

rho=1000./(1.02714+7.04e-4.*t+9,66e-7,*t.A2+3.0e-9.*t.A3-4.91e-4*(p/100)-2.633e- 

6*(p/100).*t-4.04e-9*(p/100).*t.A2-8.8e-ll*(p/100).*t./'3+1.471e-6*(p/100)A2+9.2e- 

9*  (  p/  100)  A2.  *t .  A2-3.  6  3  e  -  9  *  (  p/ 100)  A3-l.  6  6  e  -  11*  (  p/  100)  A3.  *t ) ; 

c=1570*(  1.  000-2.  15e  -  3.  *t  +4.  Oe  -  6.  *t .  A2+2.  5e  -  6*(  p/  100)  .  *t  +2.  22e  -  3*(  p/  100)  -  3.  Oe- 

6*(  p/  100) A2 ) ; 

%%  Calculation  density  and  speed  of  sound  sea  water 
[  r  h o _ s ]  =  gsw  r  ho_t_exact  ( SA,  t ,  p)  ; 

[ c  _  s ] =  gsw_sound_speed_t_exact  (  SA,  t ,  p); 

%%  r  h  o  *  c  products 

rho_c_c=rho.  *c;  %  Characteristic!  mpedance 

Castor  oil 

r  h  o  _  c  _  s  =r  h  o  _  s .  *  c_s ;  %  Char  act  er  i  st  i  c  Impedance 

Seawater 

del  ta=(  ( abs(  rho_c_c  -  rho_c_s) ) .  /  rho_c_s)*100;  %  Difference  in  Impedance  [  %] 

%%  Plot 

figure  (1) 

subplot  (2,1,1) 

pi  ot  ( t ,  r  ho) 

hoi  d  on 

grid  on 
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pi ot ( t , r  ho_s , '  r ' ) 

xl  a  bel  | '  Temper  t  ur  e  [degrees  Cel  c  i  us] '  ) 
yl  abel  ( '  Densi  ty  [kg  /  mA3]  ’ ) 

hi  e  g  1  =  I  egendf '  castor  oi  I  ' , '  sea  water' ,  '  Locati  on' ,  '  South  West'  ) ; 

t  i  1 1  e( 1  Dens  i  t  y  and  sound  speed  in  Castor  Oil  and  Seawater  at 
10m’ , '  FontWei  ght' , '  bol  d’ ) 

subplot  (2,1,2) 

pi  ot ( t ,  c) 

hoi  d  on 

grid  on 

pi  ot ( t ,  c _ s ,  1  r ' ) 

xl  a  bel  ( '  Temper  t  ur  e  [degrees  Cel  c  i  us] '  ) 
yl  abel  ('  speed  of  sound  [m/s]') 

hi  egl=l  egend( '  castor  oil', 'sea  wa  ter’, ’Location’, ’South’); 

f  i  g u r  e ( 2) 

subplot  (2,1,1) 

plot(t,rho_c_c) 

hoi  d  on 

grid  on 

pi  ot ( t , r  ho_c_s ,  '  r '  ) 

xl  a  bel  ( '  Temper  t  ur  e  [degrees  Cel  c  i  us  ]  ’ ,  ’  Font  Si  ze’ ,  16,  '  Font  Wei  ght  ’ ,  ’  demi  ) 

yl  abel  (’  I  mpedance  [  Pa*s/m] ' , '  FontSi  ze’ ,  16,  ’  FontWei  ght' , '  demi  '  ) 

hi e g 1 = I  egendf ' castor 

oi  I  '  ,  '  seawater'  ,  '  Locati  on'  ,  '  South  We  st’ ,  '  FontSi  ze’  ,  16,  '  Font  We  ight'.'demi'); 

t  i  1 1  e( '  Char  act  er  I  st  i  c  Impedances  of  Castor  Oil  and  Seawater  at 

10m'  ,  '  FontSi  ze'  ,  20,  ’  FontWei  ght'  ,  '  bol  d’  ) 

subpl  ot  ( 2,  1, 2) 

pi  ot  ( t ,  delta) 

grid  on 

xlabel(‘Temperture  [degrees  Cel  ci  us] ' ,  ’  FontSi  ze’ ,  16, '  FontWei  ght  ’ , '  demi  '  ) 
yl  abel  ('  Difference  in  Impedance  [%]’,’  Font  Si  ze' ,  16,  '  Font  Wei  ght demi  1 ) 
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APPENDIX  D.  MATLAB  CODE  FOR  PERFORMANCE 
PREDICTION 

%%  Sensitivity  prediction 

%  Script  determines  the  temperature  differences  between  the  wires  for  both  air  and 
Castor  oil  based  on  the 

%  the  formulae  in  the  article  by  J.W.  Honsbroek,  "Analytic  model  of  a  two -wire 
t  h  e  r  ma I 

%  sensor  for  flow  and  sound  measurements).”  It  also  mo  dels  the 
%  approximation  given  in  the  articles.  All  parameters  are  made  temperature 
%  dependent  to  predict  performance  over  a  wide  range  of  temperatures. 

%%  Clean  up 
cl  ear; 
cl  c; 

close  all 

f  o  r  ma  t  long 

%%  Basic  p  a  r  a  me  t  e  r  s 

f  =0:  1;  1  0  0  0; 

T  =3  1  3:  1:  6  7  3; 

TT  =4  0 :  1:4  0  0; 

TTT  =6  7  3 ; 

p  =1  0 . 1  3  2  5; 

st  ep=20e-  6; 

v _ a  =0 .  0  2  4  2  0; 
air)  [  ml  s  ] 

v _ c  =6 .  1 8 6 e -  06; 

[  m/  s  ] 

[  X  Y] =mes  hgr i  d( f ,  T) ; 

%  Mi  c  r  of  I  own 
I  y  =0 .  0  01; 

L  =2 e -  6; 

A  =2  e  -  10; 

I  0=9.  8 7 1  e -  2; 
a_a=100e-6; 


%  frequency  range  [Hz] 

%  T  e  mp  e  r  a  t  u  r  e  range  [K] 

%  T  e  mp  e  r  a  t  u  r  e  range  in  Celcius 

%  a  p  p  r  o  x  i  ma  t  e  operating  temperature  (air)  [deg  C] 

%  At  mo  spheric  Pressure  in  d  B  a  r 
%  step  distance  between  wire  [  m] 

%  a  b  i  t  r  a  r  y  particle  velocity  (IL  78  dB  re  le-12  W/  mA2  in 

%  Particle  velocity  (IL  78  dB  re  le-12  W/  mA2  in  Castor  Oil) 


%  Length  of  the  f  i  I  a  me  n  t  s  [  m] 

%  E  s  t  i  ma  t  e  d  half  width  of  the  f  i  I  a  me  n  t  [  m] 

%  Crosssection  of  a  filament 
%  Current  based  on  P  =  1 2 mW  i  n  air  [A] 

%  distance  between  wire  in  the  microflown  Titan  [  m) 
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a_c=20e-6;  %  initial  (variable)  distance  between  wire  potential 

hy dr  of  I  own  [  m] 

a_max=100e  -  6;  %  maxi  mum  di  st  ance  between  wires 

ns  =a_  max/  s  t  e  p ;  %  Number  of  steps 

h_pt  =100e-  9;  %  thickness  of  pi  at  i  num  I  ayer  [  m] 

al  pha=3.  93e  -  3;  %  Temperature  coefficient  of  resistivity  Platinum  [CA-1] 

h_Si  Ni  =200e-9;  %  thickness  of  SiNi  layer  [  m] 

r  hoc_pt  =2. 85e6;  %  density  times  specific  heat  Pt  []/(nT3*K)l 

rhoc_Si  Ni  =1. 66e6;  %  density  times  specific  heat  SiNi  ]J/(mA3*K)] 

%%  P  r  e  -  a  I  I  o  c  a  t  e  ma  t  r  i  c  e  s 

rati  o_FR=zeros(l  ength(T), I  ength(f)); 

rati  o _ DT  =z e r o s ( I  ength(T),  I  ength(f)); 

FR_a=zeros(l  ength(T),  I  ength(f)); 
predl=zeros(l  ength(T),  I  ength(f)); 
predll=zeros(  I  ength(T),  I  ength(f)); 

%%  Calculate  Power  to  acomplish  temperature 

resi  st=3.  675e  -  10*TT+9.  93e-S;  %  resistivity  of  Platinum 

R0=resi  st*l  y/A;  %  Resistance  of  the  filaments  over  temperature 

P=l  0A2.  *R0; 

%%  Calibration  report 
%%  input  p  a  r  a  me  t  e  r  s 
%  sensitivity 

suh_250Hz=122.  5;  %  sensitivity  (high  gain)  at  2  5  0  H  z  [mV/Pa] 

sul  _250Hz=0.  9732;  %  sensi  t  i  vi  ty  ( I  ow  gai  n)  at  2  5  0  Hz  [mV/Pa] 

%s  e  n  s  i  t  i  v  i  t  y  corner  frequencies 
f c 1 u  =7  7 ;  %  [  Hz] 

f  c 2 u  =7 9 1 ;  %  [  Hz] 

f  c3u=1  8  4  3  7;  %  [  Hz] 

f  c  4 u  =6 0 ;  %  [  Hz ] 

%  phase  corner  frequencies 
c 1 u  =3 0 ;  %  [  Hz] 

c  2  u  =7  0  8 ;  %  [  Hz ] 

c3u=19974;  %  [  Hz] 


c  4  u  =6 0 ; 


%  [  Hz  ] 
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%%  Calculations  velocity  e  I  e  me  n  t  sensitivity  (High  gain  /  Low  gain) 

[su  h]  =s  u  h  250Hz./((sqrt(l+(f.A2./fc3uA2))).*(sqrt(l+(fcluA2./f,A2))).*(sqrt(l+(fc4 
u A2 .  /  f .  A2)  I) .  *(  sqrt  ( l+(  f .  A2.  /f  c2uA2) ) ) ) ; 

[  su  I  ]  =s  u  I  2  5  0  Hz.  /(( sqrt  (  l+(  f  .  A2.  /  f  c3uA2)  )).*(  sqrt  (  l+(  fcluA2.  /  f  .  A2)  )).*(  sqrt  (  l+(  f  c4 
uA2.  It  .  A2)  )  )  .  *(  sqrt  (  l+(  f .  A2.  /f  c2uA2)  )  )  )  ; 

% 

%%  Calculation  of  the  properties  of  air  and  Castor  oil 
% 

%  T  h  e  r  ma  I  conductivity  air  and  Castor  oil 

k_a=-  6,  1 5 5 e -  11,  *T.  A3+6,  8 2 1  e -  8.  *T.  A2  +4.  7  4  3  e -  5.  *T-  7.  4 1 6 e -  3; 

k_ c  =5 .  5  5  6  e -  5.  *T+0.  1  6  3  7; 


%  Density  at  standard  atmospheric  pressure  [kg/mA3] 
r  h o_  a  =-  7.  4  9  9  e -  9.  *T.  A3+l.  4  5  2  e -  5.  *T.  A2  -  0.  0  1  0  2  8.  * T  +3 .  153; 

r  ho  c  =  1  0  0  0 .  /  (  1.  02714+7.  0  4  e -  4.  * T T  +  9 .  6 6 e -  7.  *  TT .  A2  +  3 .  Oe  -  9.  *  TT .  A3-4.  9 1  e -  4*  (  p  /  100)  ■ 

2.  6  3  3  e  -  6  *  (  p  / 100)  .  *TT  -  4.  0  4  e  -  9*  (  p  /  100)  .  *TT,  A2-8.  8e  -  1 1  *  (  p/  100)  .  *  T  T .  A3+l.  4  7 1  e  - 

6*(  p/  100)  A2+9.  2e  -  9*(  p/  100)  A2.  *TT.  A2-3.  6 3 e -  9*(  pi  100)  A3-l.  6 6 e -  11*  (  p/  100)  A3.  *TT)  ; 

% . 

%  Specific  Heat  of  air  (  4  0  0  deg  C)  and  Castor  oil  (40  •  4  0  0  deg  C) 

%  Formula  for  the  approximation  of  the  specific  heats 

c  p_  a  =7 .  8  9  4  e  -  11,  *T.  A5  - 1.  9 1 6  e  -  7.  *T.  A4+l.  8  2  0  e  -  4 .  *T,  A3-8.  4 1 7  e  -  2 .  *T.  A2+19.  04.  *  T  -  6  8  6.  8; 
c  p_ c  =0 .  0416.  *T.  A2  -  1  5.  3  3  6.  * T  +2  6  8 4 .  4; 

%%  Cal  cul  ati  on  of  partial  p  a  r  a  me  t  e  r  s 

D_a=k_a.  /  ( rho_a.  *cp_a) ;  %  T  h  e  r  ma  I  diffusion  coefficient  air 

D_c=k_c.  /  ( r  ho_c.  *cp_c) ;  %  T  h  e  r  ma  I  diffusion  coefficient  Castor  oil 


f  hc_a=D_a.  *rho_a.  *cp_a.  /  ( 2*pi  *L*(  200e  -  9*rhoc_Si  Ni  +100e  -  9*rhoc_pt) ) ;  %  corner 

f  requency  air  f  Hz] 

fhc_c=D_c.*rho_c.*cp_c./(2*pi*L*(200e-9*rhoc_SiNi+100e-9*rhoc_pt));  %  corner 

frequency  Castor  oi  I  [ Hz ] 


f  b a r _ L a  =(  2*pi  *LA2.  /  D_a)  ’  *f ; 
f  b  a  r  _  L  c  =(  2*pi  *LA2.  /  D_c)  ’  *f ; 
f  bar_Aa=(  2*pi  *a_aA2.  /  D_  a ) 1  *f; 
f  d_ a  =D_ a .  /  ( 2*  pi  .  * ( a_ a  A2) ) ; 


%  d i  me n s  i  o n I  e s s  frequency  at  L  (air) 

%  di  mens  i  onl  ess  frequency  at  L  (Castor  oil) 

%  d i  me n s  i  o n I  e s s  frequency  at  a  (air) 

%  Second  corner  frequency  (air) 


Calculation  for  air  with  constant  distance  between  the  filaments 
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%  Calculating  the  factor  matrix  befor  the  Bessel  function 
C_a=(((v_a*a_a.*P)./((li*pi*ly).*D_a.*k_a))'*exp(li*2*pi.*f))./fbar_Aa;  % 

% . 

%  N  u  me  r  a  t  o  r  of  the  Bessel  function  part 

numbes_a=l-(sqrt(li.*fbar_Aa).*besselk(l,sqrt(li.*fbar_Aa))); 

denombes_a=l  -  ( (  li  .  /  fhc_a’  )*f) .  *(  bessel  k(  0,  sqrtf  li  .  *fbar_La)  )- 
bessel  k(0,  sqrt(  li  .  *fbar_Aa) ) ) ; 

% . 

%  The  temperature  difference  between  the  filaments 

DT_a=C_a.  *  ( numbes_a.  /  d  e  n  o  mb  e  s  _  a ) ;  %  DT  for  filaments  in  air 

%%  Cal  cul  ati  on  of  Frequency  response  in  air 

DT  a  0  =(  (  (  v  a*  a  a .  *  P )  .  /  ( (  2*  pi  *  I  y )  .  *  D  a.*k  a)  )  )  .  *(  I  og(  (  pi  .  *a  a )  .  /  (  2*  I  y )  )  +0 .  5  7  7  )  ;  % 

DT  at  f=0  Hz 

n  =1 ; 

FR_a=zeros(  I  engthf  DT_a_0) ,  I  ength(f) ) ; 
for  n=l:  I  ength(T) 

FR  a  ( n ,  :  )  =DT  a  0(  n )  *  ( (  1 .  /  s  q  r  t  (  1  +( ( 1 .  /  ( f  he  a  ( n ).  A2  '))*(  f .  A2  )))).*(  1 ./ s  q  r  t  (  1  +( ( 1 ./( f  d 
_a[ n) .  A2 1 ) ) *[f  7 A2) ) ) ) ) ; 

n  =n  +1 ; 

end 

%%  Calculations  for  Castor  oil  at  various  distances  between  the  filaments 
FR_c=zeros(  I  engthf  DT  _  a  _  0 ) ,  I  ength(f)); 
n  =1 ; 

for  n  =1 :  ns 

f  bar  _Ac  =(  2*  pi  .  * a_c  *2 .  /  D_c )  ’  *f ;  %  di  mensi  onl  ess  frequency 

at  a  (Castor  oil) 

f d_ c =D_ c . / ( 2* p i . * ( a_ c .  *2 ) ) ;  %  Second  corner  frequency 

(air) 

C_c=(((v_c*a_c.*P)./((li*pi*ly).*D_c.*k_c))’*exp(li*2*pi.*f)),/fbar_Ac; 

n  u  mb  es_c=l  -  ( sqrt  (  li  .  *f  bar_Ac) .  *bessel  k(  1,  sqrt{  li  .  *fbar_Ac) ) ) ; 

denombes_c=l  -  ( ( li  .  /  fhc_c'  )*f) .  *(  bessel  k(  0,  sqrtf  li  .  *fbar_Lc)  )- 
bes s el  k(  0,  s qr  t  ( li  .  * f  bar  _Ac I) ) ; 

DT_  c  =C_  c .  *  ( numbes_c.  /  d  e  n  o  mb  e  s  _  c )  ;  %  DT  for  filaments  in 

Castor  oil 

DT_c_0=(  ( (  v_c*a_c.  *P) .  /  ( ( 2*pi  *1  y) .  *D_c.  *k_c) ) ) .  *(  I  og(  ( pi  .  *a_c)  .  /  ( 2*1  y) )  +0.  577) ; 
%%  Calculate  the  frequency  response  in  Castor  oil 
n  n  =1 ; 
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for  n  n  =  1 :  I  engt  h( DT  _  c  _  0 ) 


FR  c  (  n  n ,  :  )  =DT  c  0(nn)*((l./sqrt(l+((l./(fhc  c  (  n  n ).  A2  '))*(  f .  A2  )))).*(  1 ./ s  q  r  t  (  1  +(  (  1 .  / 
(fd.c(nn).  A2'))*(f .  A2)))J); 

n  n  =n  +1 ; 

end 

n  n  =1 ; 

%  Loop  to  compare  performance  at  temperatures  to  performance  in  air 

for  n  n  =  1 :  I  engt  h( DT  _  a  _  0 ) 

rati  o_DT(nn,  :  )  =DT_  a  (  3  61,  :  )  ,  /  DT  _  c  (  n  n ,  :  ) ; 

n  n  =n  +1 ; 

end 

%%  Calculation  of  the  ratio  of  the  frequency  response 
n  n  =1 ; 

%  Loop  to  compare  performance  at  temperatures  to  performance  in  air 

for  nn=l:  I  ength(T)  %  Loop  to  compare  performance  at  temperatures  to  performance 
i  n  ai  r 

rati  o_FR(nn,  :  )=FR_a(361,  :  )./FR_c(nn,  :  ); 
rati  o  _  F  R  _  c  o  I  d(  nn,  :  )  =FR_a(  274,  :  )./FR_c(  nn,  :  ); 
n  n  =n  n  +1 ; 
end 
n  n  =1 ; 

for  n  n  =  1 : I engt  h( DT  _  a  _  0 ) 

predlf  nn,  :  )  =su_h.  /  abs(  rati  o_FR(  nn,  :  ));  %  predicted  performance  at  313  deg  C 
using  calibration  report 

predllfnn,  :  )=su_h./abs(rati  o_DT(nn,  :  ));  %  predicted  performance  at  313  deg  C 
using  calibration  report  and  raw  DT 

n  n  =n  n  +1 ; 

end 

%%  Get  e  x  p  e  r  i  me  n  t  a  I  results 

d  a  t  a  =x  I  s  r  e  a  d  ( ‘  C:\Users\Marni  x  \  Document  s  \  MATLAB  \  Thes  i  s\Fi  nal 
versi  o  n  s  \  Graphs!  Predi cti  ons\test_2. xl s  ) ; 

S  v  =[  dat  a( :  ,  1)  d  a  t  a  ( :  ,  2 )  ] ;  %  get  data  on  velocity  sensor  only 

S  v  (  :  ,  2)  =10.  A(  ( Sv( :  ,  2) .  /  20)  +9) ;  %  Response  in  [mV/Pa] 

% . 

%  Plot  prediction  at  different  distances  between  the  wires 
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f  i  g u  r  e (  100+n) 

plot(f,predl(274,:))  %  313  degrees  C 

hoi  d  on 

pi  ot(f,  predl(  161, :  ),  '  r  )  %  2  0  0  degrees  C 

s  cat  t  er ( Sv( : ,  1) ,  Sv( :  ,  2) . / 10) ; 

xl  abel  ('  Frequency  [  Hz] ' ,  ’  FontSi  ze' ,  14) 

yl  abel  ('  Sensi  ti  vi  ty  [  mV  /Pa] '  ,  ’  FontSi  ze’  ,  14) 

set(fi  ndobj  (  type', ’axes'), ’fontsize',  12, ’Font  We  ight'.'de  mi’) 

a x i  s ( [  0  1  0  0  0  0  0.5  ]) 

I  egend('  313  deg  C','200  deg  C' , ' experi  ment'  ,  '  313  deg  C  100 
mi  crons' , '  Locat  i  on' ,  '  Nort  h E a s t '  ) 

tt  =  strcat]  {‘  Predi  cted  sensitivity,  no  gain,  20%  loss  for  },  num2str(a_c*le6),{ 
mi  crons'  } ) ; 

ti  tl  e(tt,  '  FontSi  ze'  ,  18,  '  Font  Wei  g  h  t '  ,  '  bol  d  ) 
f  i  gure(  200+n) 

plot(f,predll(274,:))  %  313  degrees  C 

hoi  d  on 

pi ot(f,  predlli  161, : ), '  r' )  %  200  degrees  C 

s  cat  t  er  (  Sv(  :  ,  1) ,  Sv(  :  ,  2) .  /  10) ; 

xl  abel  ('  Frequency  [Hz]', 'FontSize' ,14) 

yl  abel  ('  Sensi  ti  vi  ty  [  mV  /  Pa] ' ,  '  Font  Si  ze’  ,  14) 

set(fi  ndobj  (  type', ’axes’), 'fontsize' ,12,'Font  We  ight'.'de  mi’) 

a x i s ( [ 0  1  0  0  0  0  5  ] ) 

I  egend]'  313  deg  C,'200  deg  C' ,  '  experi  ment’ ,  ’  313  deg  C  100 
mi  crons' ,  '  fontsi  ze' ,  12,  '  Locati  on' , '  Nort  hEast '  ) 

tt  =  strcat]  {'  Predi  cted  sensitivity  1 4  0  d  B  using  DT,  no  gain,  20%  loss  for  }, 
n  u  m2  s  t  r  ( a  _  c  *  1  e  6 ) ,  {  mi  c  r  o  n  s  } ) ; 

t  i  1 1  e  ( 1 1 ,  'FontSize',  18,  '  Font  We  i  g  h  t ' ,  'bold  ) 

% . 

figure  (1)  %  Plot  ratio  of  DT  for  various  distances  between  the 

wi  r  es 

pi ot(f, predi]  274, :  )) 
a x i s ( [ 0  1  0  0  0  0  0.5  ]) 
hoi  d  al  I 

xl  abel  ('  Frequency  [Hz]',  'Fontsize',  16) 
yl  abel  ('  Sensi  ti  vi  ty  [  mV  /Pa]',  '  Fontsize'  ,  16) 
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I  e  g  e  n  d  ( 1  2  0  u  m'  ,  '  40um’  ,  '  6  0  u  m'  ,  '  80  um'  ,  '  lOOum'  ,  '  fontsi  z  e '  ,  14,  ’  Locati  on’  ,  '  East  Outsi  d  e '  ) 

titlef'Predicted  sensitivity,  no  gain  at  313  deg 
C ' ,  '  FontSi  z  e ' ,  20,  'Font  We  i  g  h  t 1 ,  '  bold') 

figure  (2) 

pi  ot ( f ,  pr  e d  1  ( 161, : ) ) 
a x i  s ( [  0  1  0  0  0  0  0.5  ]) 
hoi  d  al  I 

xl  abel  ('  Frequency  [  Hz] ' , '  FontSi  ze' ,  14) 
yl  abel  ('  Sensi  ti  vi  ty  [mV/Pa] ' ,  ’  FontSi  ze'  ,  14) 

I  egend(‘  20um'  , '  4  0  u  m' ,  '  60um' ,  ’  80um' , '  10  0um'  , '  Locati  on' , '  EastOutsi  d  e  ’  ) ; 

t  i  1 1  e( '  Pr  edi  ct  ed  sensitivity,  no  gain  at  2  0  0  d  e  g 
C'  ,  '  FontSi  ze'  ,  20,  Font  Wei  g  h  t  ’  ,  '  bol  d  ’  ) 

figure  (3)  %  Plot  ratio  of  DT  for  various  distances  between  the 

wi  r  es 

pi  ot  ( f ,  predll]  2  74,  :  ) ) 
a x i  s ( [  0  1  0  0  0  0  0.3  ]) 
hoi  d  al  I 

xl  abel  ('  Frequency  [  Hz] ' ,  '  FontSi  ze' ,  16) 
yl  abel  ('  Sensi  ti  vi  ty  [mV/Pa] '  ,  '  FontSi  ze'  ,  16) 

I  egendf'  2  0  u  m'  ,  '  4  0  u  m’  ,  '  6  0  u  m’  ,  ’  8  0  u  m'  ,  '  100um'  ,  '  fontsi  ze'  ,  14,  '  Locati  on1  ,  1  EastOutsi  d  e '  ) 

t  i  1 1  e( '  Pr  edi  ct  ed  sensitivity  using  DT,  no  gain  at  313  deg 

C’, '  FontSi  ze',  20,  'Font  We  ight’,’  bold') 

figure  (4) 

pi  ot  ( f ,  pr  e d 1 1  ( 161, : ) ) 
a x i  s ( [  0  1  0  0  0  0  0.3  ]) 
hoi  d  al  I 

xl  abel  ('  Frequency  [  Hz] ' , '  FontSi  ze' ,  14) 
yl  abel  ('  Sensi  ti  vi  ty  [mV/Pa] '  ,  ’  FontSi  ze'  ,  14) 

I  egend] '  20um'  ,  '  40um'  ,  '  60um’  ,  '  80um'  ,  '  lOOum'  ,  '  Locat  i  on'  ,  '  EastOutsi  de'  ) ; 

ti  ti  e(  ‘  Predi  cted  sensitivity  using  DT,  no  gain  at  2  0  0  d  e  g 
C',  '  FontSi  ze',  20, 'Font  We  ight', 'bold') 

% 

figure  (300  +n) 

i  ma  g  e  s  c  ( f ,  T  (  1 :  2  7  4  )  -  2  7  3 ,  p  r  e  d  1  (  1 :  2  7  4 ,  :  )  ) 
xlabel  (‘  Frequency  [  Hz]  ’  ,  '  FontSi  ze'  ,  14) 
yl  abel  ('Temperature  [deg  C] '  ,  '  FontSi  ze'  ,  14) 
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colorbar('vert’) 
c  a x i  s  ( [  0  .  3] ) 
h  =c o I  orbar; 

yl  abel  ( h,  ’  mV/  Pa’ , '  FontSi  z e ’  ,  14) 

tt  =  s  t  r  c  a  t  ( { '  Sensitivity  prediction  vs  filament  temperature  for  }, 
n  u  m2  s  t  r  ( a  _  c  *  1  e  6 ) ,  {  mi  c  r  o  n  s  } ) ; 

ti  1 1  e(tt,  ’  FontSi  ze' ,  18,  '  FontWei  ght' , '  bol  d' ) 

figure  (400  +n) 

i  ma  g  e  s  c  ( f ,  T  (  1 :  2  7  4  )  ■  2  7  3 ,  p  r  e  d  1  1  (  1 :  2  7  4 ,  :  ) ) 
xl  abel  ( '  Frequency  [  Hz] ' , '  FontSi  ze' ,  14) 
yl  abel  ('  Temperature  [deg  C] '  ,  ’  FontSi  ze’  ,  14) 
colorbar('vert') 
c a x i  s ( [  0  . 3] ) 
h  =c o I  orbar; 

yl  abel  ( h,  '  mV/  Pa’  ,  '  FontSi  ze’  ,  14) 

tt  =  strcat(  {'  Sensi  ti  vi  ty  prediction  vs  filament  temperature  using  DT  for  }, 
num2str(a_c*le6),{  mi  crons  }); 

ti  ti  e(tt,  ’  FontSi  ze'  ,  18,  ’  Font  We  i  ght' ,  '  bol  d' ) 


n  =n  +1 ; 

a_c=a_c+st  ep; 
end 
n  =1 ; 

for  n  =1 :  ns 
fi  g  u  r  e  (  100+n) 
hoi  d  on 

pi  ot  ( f ,  pr  edl(  2  7  4, ;  )  ,  '  -  b '  ) 

I  egend( ' 313  deg  C','200  deg  C’ ,  ’  experi  ment' , ' 313  deg  C  100 
mi  crons’  ,  ’  Locati  on'  ,  '  NorthEast'  ) 

n  =n  +1 ; 

end 

for  n  =1 :  ns 
f  i  gur  e(  2  0  0  +n) 
hoi  d  on 


pi  ot  ( f ,  pr  edll(  2  7  4,  :),’••  b’ ) 
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I  egend('  313  deg  C '  ,  '  2  0  0  deg  C’  ,  '  experi  merit’ ,  ’  313  deg  C  100 
mi  crons' ,  ’  Locat  i  on’ ,  '  Nort  h E a s t '  ) 

n  =n  +1 ; 

end 

%%  predi  cti  ons  at  313  deg  Cel c i  us 
f  i  g u r e( 5) 

pi  ot ( f ,  10.  *  pr  ed  1  (  2  7  4, :  ) ) 
hoi  d  on 

pi  ot  (  10.  * p r  edll(  2  7  4,  :  ) ,  '  r '  ) 
s  cat  t  er  ( Sv( : , 1 ) , Sv ( : ,  2 ) ) ; 

legend!' prediction  using  approximation',1  predi  cti  on  using  DT  ratio','  averaged 
experi  mental  resul  t  s ' , '  Locati  on'  ,  '  Nort  hEast ’ ) 

xl  abel  ('  Frequency  [  Hz] '  ,  '  FontSi  ze'  ,  16) 

yl  abel  ('  Sensi  ti  vi  ty  [  mV  /  Pa] ' ,  '  Font  Si  ze’  ,  16) 

set(fi  ndobj  (  type’  ,  '  axes'  ),  'fontsi  ze'  ,  14,  '  FontWei  ght'  ,  '  demi  '  ) 

a x i s ( [ 0  1  0  0  0  0  1  ]) 

ti  tl  e('  Sensi  ti  vi  ty  1 4  0  d  B ,  no  gain,  20%  loss'  ,  ’FontSi  ze',  20, 'Font  We  ight'.'bold') 
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